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I. 


INTRODUCTION 


There  is  considerable  interest  in  the  development  of  extremely 
sensitive  superconducting  magnetometers  for  the  detection  of  changes 
in  magnetic  field  intensity.  The  most  sensitive  detectors  developed 
up  to  this  time  are  thin  film  dc  SQUID  magnetometers!  The  dc  SQUID  con- 

9 

sists  of  a  superconducting  loop  interrupted  by  two  Josephson  junctions. 

If  the  junctions  are  biased  near  the  critical  current  by  an  externally 
applied  dc  bias  current,  the  voltage  across  the  junctions  becomes  a  sen¬ 
sitive  function  of  the  flux  threading  the  loop.  Thus,  the  dc  SQUID  acts 
like  a  flux  to  voltage  transducer. 

The  resolution  of  the  dc  SQUID  magnetometer  is  determined  pre- 
dominately  by  three  factors.  First,  the  Josephson  junctions  employed 
in  the  device  must  be  non-hysteretic.  As  a  result,  the  thin  film  tunnel 
junctions  typically  used  in  the  device  are  externally  shunted  by  a  normal 
resistance.  The  Johnson  noise  associated  with  this  resistance  produces 
fluctuations  in  the  voltage  across  the  SQUID.  ,  This  noise  can  be  referred 
to  the  input  as  an  effective  flux  noise,  Sc.  A  useful  figure  of  merit 

is  then  the  energy  factor  =  Sj/2L,  where  L  is  the  SQUID  loop  inductance. 

In  the  usual  mode  of  operation,  the  incident  magnetic  field  is  coupled 
into  the  SQUID  through  a  superconducting  transformer.  The  pick-up  coil 
is  connected  in  a  series  with  a  coupling  coil  which  is  inductively  coupled 

to  the  SQUID  loop.  In  this  case,  the  minimun  detectable  change  in  magnetic 

field  energy  stored  in  the  pick-up  coil  is  proportional  to  the  energy 
factor  S^.  For  conventional  dc  SQUIDs,  the  energy  factor  is  determined 
by  the  temperature  T,  the  loop  inductance  (_,  and  junction  capacitance  C, 

Sr  a  kDT  AX.  The  ultimate  resolution  for  quantum  noise  limited  devices 
appears  to  be  S^  %  h.b 

The  second  factor  which  limits  the  resolution  of  dc  SQUID  magneto- 

3 

meters  is  the  presence  of  1/f  noise  at  the  output  of  the  device.  The 
source  of  this  noise  is  not  well  understood.  However,  it  is  likely  that 
the  dominant  source  is  related  to  low  frequency  thermal  fluctuations 
within  the  Josephson  junction.^’0  These  fluctuations  increase  in  importance 


as  the  effective  junction  volume  is  reduced.  Since  the  effective  volume 

scales  with  the  junction  capacitance,  attempts  to  reduce  the  energy  factor 

at  low  frequencies  by  decreasing  the  junction  capacitance  have  failed. 

In  fact,  for  the  most  sensitive  devices  with  Sr  ^  h,  1/f  noise  is  the 

>  £  o 

dominant  noise  source  at  frequencies  f  'v  10-100  KHz. 

The  third  factor  which  limits  the  resolution  of  the  dc  SQUID  magneto¬ 
meter  is  the  intrinsic  low  loop  inductance  of  the  device.  For  example, 
the  maximum  practical  loop  inductance  for  devices  operated  at  4K  is  L  ^ 

nH.  SQUIDs  with  energy  factors  Sr  approaching  h  have  been  fabricated  with  loop 
9  L 

inductance  L  ^pH.  However,  the  magnetic  field  resolution  depends  not 
only  on  the  energy  factor  SE,  but  also  on  the  effective  area  of  the  pick¬ 
up  coil.  Since  the  pick-up  coil  inductance  must  be  matched  to  the  induc¬ 
tance  of  the  coupling  coil  to  the  SQUID,  the  magnetic  field  resolution 
is  limited  by  the  coupling  efficiency  between  the  SQUID  loop  and  the 
coupling  coil.  As  a  result,  planar  dc  SQUIDs  with  loop  inductance  L' 
pH  show  little  promise  as  sensitive  detectors  of  magnetic  field  intensity. 

Recently,  a  novel  superconducting  device  related  to  the  dc  SQUID 
has  been  suggested  by  Ho  Jung  Paik.10  In  this  device,  the  single  dc  SQUID 
loop  is  replaced  by  a  double  loop  structure.  In  this  configuration, 
a  "SQUID"  loop  containing  two  Josephson  junctions  and  a  superconducting 
"coupling"  loop  are  connected  in  parallel  with  a  shorting  capacitance. 

It  has  been  suggested  that  the  behavior  of  this  device  will  be  closely 
related  to  the  ordinary  dc  SQUID.  In  particular,  a  heuristic  calculation 
suggests  that  the  coupling  properties  of  the  device  will  be  determined 
by  the  coupling  loop  inductance,  whereas  the  energy  factor  will  be  determined 
by  the  SQUID  loop  inductance.  Thus,  by  decreasing  the  SOUID  loop 
inductance  without  sacrificing  the  ability  to  couple  efficiently  to 
larqe  input  coils. 

In  this  report  we  describe  the  results  of  a  detailed  investigation 
of  the  double  loop  SQUID  magnetometer.  We  begin  by  presenting  a  thermal 
activation  model  for  the  DC  SQUID  which  supplements  the  original  numerical 
simulations  Of  the  device.  The  model  is  particularly  useful  for  deriving 
the  behavior  of  the  device  at  low  bias  currents.  The  equations  of  motion 


for  the  dc  SQUID  are  then  compared  with  a  lumped  circuit  element  model  for 
the  double  loop  SQUID.  A  relationship  between  the  double  loop  and  con¬ 
ventional  SQUID  characteristics  in  the  absence  of  noise  is  described.  Direct 
numerical  simulations  of  the  double  loop  device  are  then  performed.  Noise 
free  characteristics  are  computed  as  a  function  of  the  coupling  loop  resonant 
frequency,  the  bias  current  and  the  applied  flux.  Noise  rounded  character¬ 
istics  are  then  generated.  An  approximate  expression  for  the  flux  resolution 
as  a  function  of  the  total  loop  inductance  is  determined.  A  model  is 
developed  for  conventional  and  double  loop  SQUID  magnetometers.  The  model 
includes  explicitly  the  effects  of  the  input  circuit  on  the  SQUID  character¬ 
istics.  The  resolution  of  superconducting  and  lossy  double  loop  magnetometers 
is  calculated  as  a  function  of  the  signal  frequency  and  the  input  and  SQUID 
circuit  parameters. 

Some  of  the  basic  results  reported  for  the  isolated  double  loop  SOUID 
were  computed  in  conjunction  with  a  project  on  the  analysis  of  double  loop 
linear  amplifier  performed  for  the  National  Science  Foundation  under  Grant 
DAR-8009388. 


II.  MODEL  FOR  THE  CONVENTIONAL  DC  SQUID 

A.  Numerical  Simulation  of  Lumped  Circuit  Element  Model 

The  conventional  dc  SQUID^consi sts  of  a  superconducting  loop  inter- 

2 

rupted  by  two  Josephson  junctions.  A  lumped  circuit  element  model  for 

the  device  is  shown  in  Figure  1?  The  SQUID  is  biased  at  constant  current 

I.  The  average  voltage  V  across  the  SQUID  output  is  then  a  periodic 

function  of  the  flux  !’  applied  to  the  loop  inductance  L.  The  RSJ  model  (Fig.  2(a)) 

a  ^2 

is  used  to  describe  the  behavior  of  the  junctions.  The  junction  critical 
currents  are  IQl  =  (1  -  a)IQ  and  Iq2  =  (1  +  a) IQ ,  where  a  is  the  asymmetry 
parameter.  The  junction  shunt  resistance  is  R.  We  consider  the  limit 
Bc  =  RC/(4>0/2ttIoR)  «  1. 

The  differential  equations  describing  the  development  of  the  quantum 

mechanical  phase  differences  6,  and  across  the  Josephson  junctions 
4  1  t- 

are  then 


^  d6 


H'  UV.U 

ar ■>«  -  j  - 10  a  - «)  «,  ♦  vw1/r. 


d6? 

ST  =  1/2  +  J  "  I0{1  +  a)  Sln  62  +  VN2/R’  (2) 

where  and  are  the  Johnson  noise  voltages  associates  with  the 
shunt  resistances.  The  circulating  current  J  is  related  to  the  junction 
phase  differences  by  the  quantization  of  the  fluxoid 

LJ  =  (.Sj  -  ,S2)/2n  -  -i'fl  (3) 

where  <J>  is  the  quasistatic  flux  applied  to  the  ring.  The  voltage  developed 
acorss  the  junctions  is 

4> 

V 


t 

f 
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The  equations  of  motion  have  been  solved  by  nimerical  simulation?’ 

The  voltage  V(t)  and  circulating  current  J(t)  are  computed  as  a  function 

of  the  bias  current  I  and  applied  flux  >i>  by  stepwise  integration  of  Eqs. 

1-4.  The  time  averaged  voltage  V(l, ),  forward  transfer  function  'sV/at  , 

a  a 

voltage  noise  power  spectral  density  Sv(f),  circulating  current  noise 
spectral  density  Sj(f),  and  correlation  noise  spectral  density  Syj(f)  are 
then  computed.  The  effective  flux  noise  spectral  density  at  the  input  is  S 
Sw/OV/W  ).  In  most  cases,  we  are  interested  in  the  constant  low  fre¬ 
quency  component,  S°,  at  frequencies  f  <<  IQR/<;  . 

The  numerical  simulation  yields  results  which  are  in  good  agreement 
with  observed  characteristics  for  bias  currents  I  o  21  .  However,  the 
simulation  is  fairly  tedious.  In  addition,  at  bias  currents  I  <<  21, 
numerical  simulation  is  not  practical.  Furthermore,  the  simulation 
does  not  provide  a  clear  intuitive  framework  upon  which  further  investiga¬ 
tion  can  be  based. 

B .  Thermal  Activation  Model 

We  present  a  thermal  activation  model  for  the  conventional  dc  SQUID 
which  provides  a  clear,  intuitive  description  of  the  behavior  of  the  device 
The  model  depends  on  the  fact  that  the  equations  of  motion  (Eqs.  1-4) 
can  be  derived  from  a  potential  U{6^,6,,).  This  is  not  a  trivial  statement 
In  fact,  the  equations  of  motion  of  numerous  interesting  devices  cannot 
be  related  to  a  potential  field.  The  double  loop  SQUID  is  an 
example  of  such  a  device.  However,  both  the  single 
resistively  shunted  Josephson  junction  and  the  rf  SQUID  can  be 
described  by  a  potential^’ ^ For  example,  the  equations  of 
motion  for  the  single  current  biased  Josephson  junction  (RSJ  model) 
are 


B  £«  +  _ 

c  d?  dr) 


d  (S 


DU 

DtS 


=  -  —■  +  V. 


(5) 


where  u  (A)  =  -  cos  A  -  i 1 . 


(6) 
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In  these  equations,  the  bias  current  I  =  i I  ,  and  the  voltage  across  the 
junction  V  =  (I0R)v  =  (IoR)d5/de.  The  dimensionless  time  0  =  t/(oo/2rIQR). 
At  bias  currents  i  <  1,  the  potential  consists  of  a  one  dimensional  chain 
of  wells  separated  by  barriers  (the  "washboard"  potential,  Fiq.  2(b)) 

develops  in  time  like  the  coordinate  of  a  particle  moving  in  a  viscous 
medium  within  the  wells.  The  random  Johnson  noise  voltages  v^  cause 
the  particle  to  make  transitions  over  the  barrier  between  the  wells. 

A  voltage  pulse  is  generated  in  the  process.  Thus,  the  time  averaged 
voltage  is  proportional  to  the  average  rate  of  escape  over  the  barriers. 

At  low  temperatures  (k^T  «  2eIQR),  the  power  spectral  density  of 
the  thermal  fluctuations  in  the  voltage  across  the  shunt  resistor  will 
roll  off  like  S^(m)  ^  ftui  exp  (-Tiw/kgT).  In  this  regime,  the  dominant 
source  of  voltage  fluctuations  across  the  single  junction  will  be  produced 
by  the  detailed  quantum  mechanical  behavior  of  the  junction.  The  inves¬ 
tigation  of  this  behavior  is  a  subject  of  consideration  interest  at  this 
time. 


The  conventional  dc  SQUID  equations  can  be  related  to  the  equations 
of  motion  of  a  particle  moving  in  a  two  dimensional  potential  field. 


"*  r 

:>c  T 

dt 


dt 


=  -  +  v  k  =  1  2 

3<S  ^  K  1  ’  c 


(7) 


where  the  potential  u(6p  <$2)  satisfies 

u(«Sj ,A2)  =  -  cos  'Sj  -  cos  <5 2  -  (^  +  $2)  i/2  +  ^Kj2/2  (8) 

and 

j  =  (<5j  -  <$2  •  27T<J>a)/TT3.  (9) 

The  junction  parameters  are  critical  current  I0,  capacitance  C  and  shunt 

resistance  R.  the  voltage  developed  across  the  junctions,  V  =  vIqR,  is 

a  function  of  the  dc  bias  current  I  =  i  I  and  the  applied  flux  4'  =  t>  4- 

0  a  a  0 


6 


The  screening  factor  is  tf  =  2 L IQ /'TQ  and  b  =  2nl  R^C/*  .  We  shall  be 
interested  in  the  solution  to  these  equations  in  the  overdamped  limit, 

*c  «  U 

Equations  7-9  are  a  set  of  classical  equations  of  motion  for  a  particle 
moving  in  a  two  dimensional  potential  field.  The  form  of  the  potential 
is  similar  to  the  one  dimensional  field  for  the  single  current 
biased  junction.  The  potential  consists  of  a  family  of  wells  located  within 
a  low  potential  trough.  An  example  is  plotted  in  Figure3  for  bias  current 
i  =  0,  applied  flux  .f»  =  0.5  and  screening  factor  0  =  1.  At  low  bias  currents 
in  the  absence  of  noise,  the  particle  velocity  will  be  damped  out.  Thus,  the 
particle  will  settle  into  a  single  well.  As  the  bias  current  is  increased, 

the  potential  trough  is  tipped  and  the  barrier  between  successive  wells  is 

decreased.  At  some  critical  current,  ic,  which  is  a  function  of  the  applied 
flux  :a,  the  barrier  vanishes.  At  this  point,  the  particle  propagates 
freely  down  the  potential  trough.  Since  the  voltage  is  proportional  to  the 
direct  sum  of  the  components  of  the  particle  velocity,  a  voltage  develops 
across  the  SQUID  (see  Fig.  4  }. 

In  the  presence  of  noise  ,  v^  f  0),  the  particle  undergoes  Brownian 
motion  in  the  potential  field.  At  bias  currents  i  <  i  (c  ),  the  particle 

motion  occurs  predomi nately  within  the  potential  well.  However,  the  particle 

will  occasionally  make  a  classical  transition  over  the  barrier.  The  thermal 
activation  of  the  particle  from  one  well  to  the  next  produces  an  average 
voltage  across  the  SQUID  even  at  bias  currents  i  <  i  (<■  )  (see  Fig.  5  ). 

The  voltage  pulse  produced  as  the  particle  escapes  from  the  well  is  usually 
accompanied  by  a  pulse  in  the  circulating  current  around  the  SQUID  as 
each  of  the  phase  differences  advances  by  2n.  Fluctuations  in  the  number 
of  pulses  per  unit  time  generate  low  frequency  noise  in  the  voltage  and 
circulating  current.  In  addition,  noise  is  generated  by  the  random  motion 
of  the  particle  in  the  wells.  This  source  of  noise  is  most  pronounced  in 
the  low  frequency  current  noise  at  bias  currents  i  <■.  i  (<j>  ).  Motion 
within  the  well  does  not  contribute  substantially  to  the  low  frequency 
voltage  noise.  An  example  of  current  and  voltage  low  frequency  noise 
spectral  densities  as  a  function  of  bias  current  are  plotted  in  Fig.  6 
and  Fig.  7.  Near  the  critical  current,  the  random  generation  of  the 
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current  and  voltage  pulses  is  the  dominant  source  of  noise.  Since  the 
current  and  voltage  pulses  are  coincident,  the  two  noise  components  are 
strongly  correlated  (see  Fig. 8  ). 

The  conventional  dc  SQUID  I-V  characteristics  and  noise  spectral 
densities  can  be  computed  directly  from  the  form  of  the  potential  u 
for  values  of  bias  current  i  ic.  The  probability  of  making  a  single 
transition  over  a  barrier  height  Au  in  the  overdamped  limit  (8C  <-•  1) 
depends  exponentially  on  the  barrier  height. 


CwlCw2Csl 


exp 


(10) 


where  Cwl’  Cw2  ’  and  Csj,  j !  are  the  curvatures  of  the  well  and  saddle 
expressed  in  normal  coordinates.  The  factor  is  the  ratio  of  the  thermal 
energy  to  the  intrinsic  energy  of  the  junction,  r  =  2rkgT/I  <?Q. 

The  motion  of  the  particle  down  the  chain  of  well  is  analogous  to  the 
action  of  a  non-paralyzable  counter.  The  occurence  of  an  event  at  the 
counter  corresponds  to  the  escape  of  the  particle  over  the  barrier.  The  counter 
deal  time  corresponds  to  the  transit  time  t  over  the  barrier.  The  times  be¬ 
tween  successive  pulses  are  the  sum  of  the  independent,  randomly  distributed 
transition  and  pulse  times.  The  transition  times  (time  in  the  wells)  are 

•  p 

exponentially  distributed  with  mean  1/p  and  variance  1/p  .  We  assume  that 

2 

the  variance  in  the  pulse  time  i  •  1/p  .  The  mean  number  of  pulses 

v  2  2  3 

arriving  in  unit  time  is  n  =  1/p  and  the  variance  in  n  is  (An)  =  /p  , 

9  . 

where  o  is  the  sum  of  the  variances  in  the  transition  time  and  p  is  the 

pp 

sum  of  the  means . 

Thus,  the  effective  transition  rate  is 


r  =  p/(l  +  pi) 


(11) 


and  the  variance  in  the  number  of  transitions  per  unit  time  is 


(An) 


2  , _ EL 


(1  +  pt  )’ 


(12) 


The  average  voltage  is 


v  =  2vr 


2np 

( 1  +  P  i ) 


j 

i  . 
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(13) 


J 


and  the  low  frequency  voltage  spectral  density  generated  by  fluctuations 
in  the  number  of  pulses  is 


S°  =  2(2tt)2  An2  =  2 ( 2tt ) 2  - 2 -  ,  f«  f  (14) 

V  (1  +  Pt)3  ° 


As  mentioned  previously,  the  circulating  current  noise,  S° ,  arises 
from  two  sources.  There  is  a  contribution  generated  by  the  motion  of  the 
particle  from  one  well  to  another,  and  a  contribution  from  motion  within 
the  well.  In  the  first  case,  the  current  pulse  jp  =  J*  (j  -  Jw)dO,  where 
the  limits  of  integration  are  taken  from  the  bottom  of  the  following  well, 
and  j  is  the  circulating  current  in  the  well.  For  average  pulse  size 
<j  ,  the  contribution  from  the  fluctuation  in  the  number  of  pulses  per 
unit  time  is 


?(y2p 

(1  +  Pi)3 


(15) 


The  contribution  to  S°  from  the  motion  within  the  well  can  be 

J 

estimated  in  the  low  noise  limit  by  ordinary  first  order  correlation 

theory  following  the  treatment  of  Likharev  and  Semenov  for  the  single 
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junction  at  i  <  ic.  For  a  particle  confined  within  the  well,  the  low 
frequency  current  noise  spectral  density  is 


S°  =  4r(a2  +  b2) 
w  (a  +  b  +  iigab)2 


(16) 


where 


a  =  (i  -  a)cos  :  mip  (17) 

b  =  (1  +  «)cos2  min  .  (18) 


We  assume  that  the  motion  within  the  well  is  essentially  uncorrelated 
with  the  noise  generated  by  the  random  production  of  the  current  pulses. 


t 

\ 


Thus,  for  effective  pulse  rate  r  and  dead  time  r. 


(1  -  ri 


(19) 


The  voltage  and  current  pulses  are  produced  simultaneously,  and 
thus  are  strongly  correlated.  Since  the  voltage  noise  is  dominated  by 
the  voltage  pulse  noise,  the  low  frequency  correlation  spectral  density 
can  be  approximated  by  S°.  ~  S°  . 

VJ  V  jp 

Note  that  Eqs.  11-19  are  only  valid  for  bias  currents  i  £  1.  At 
bias  current  i  k  1  such  that  Au  %  T,  the  approximations  used  to  derive 
Eq.  10  begin  to  break  down.  The  model  seems  to  work  fairly  well  for  values 
of  Au  %  v/2.  At  bias  currents  such  that  Au  <  T/2,  the  computed  voltage 
and  noise  spectral  densities  decrease  rapidly  and  the  model  fails.  Also 
note  that  Eqs.  11-19  have  been  derived  considering  motion  only  in  the  direc¬ 
tion  of  increasing  6.  However,  at  bias  currents  i  s  0,  the  potential 
difference  between  successive  wells  approaches  zero.  In  this  case,  the 
particle  can  make  a  transition  either  to  the  right  with  probability  pr  =  p 
generating  a  voltage  pulse  +  2n,  or  to  the  left  with  probability  pp  = 
p  exp  (-27ti /r )  generating  a  voltage  pulse  -2tt.  At  i  -  0,  the  rates  are 
equal,  and  v  =  0.  Near  i  -  0,  v  is  a  function  of  both  pr  and  pp.  A 
detailed  discussed  of  this  effect  can  be  found  in  Appendix  I. 
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III.  MODEL  FOR  THE  DOUBLE  LOOP  SQUID  MAGNETOMETER 
A.  Lumped  Circuit  Element  Model 

The  lumped  circuit  element  model  for  the  double  loop  SQUID  is  shown 
in  Figure  9.  The  superconducting  outer  loop  is  divided  into  two  parts  by 
a  shorting  capacitance  C$  with  loss  R$.  The  lower  "SQUID"  loop  of  induc¬ 
tance  L  contains  the  Josephson  junctions.  The  upper  "coupling"  loop  in¬ 
ductance  is  L$.  We  assume  that  the  current  feed  is  symmetric  and  the 
mutual  inductance  between  the  loops  is  negligible.  The  junction  parameters 
are  critical  current  IQ,  shunt  resistance  R  and  junction  capacitance  C. 

If  the  device  is  biased  at  constant  current  I,  a  voltage  V  appears  across 
the  output.  The  time  averaged  voltage  V  is  then  a  periodic  function  of 
the  flux  <t>  threading  the  superconducting  loop. 

The  equations  describing  the  time  development  of  the  junction  phase 
drops,  and  <?2>  are  identical  to  those  for  the  conventional  dc  SQUID. 

Vj  =  (<?0/2ir)(dvS1/dt)  =  R(I/2  -  J  -  I0  sin  <S?)  +  ,  (20) 


and 


V2  =  (4>0/2n)(d*2/dt)  =  R(  1/2  +  J-IQ  sin  S2)  +  VN?.  (21) 

where  J  is  the  current  circulating  in  the  "SQUID"  loop,  and  and  V^2 
are  random  voltage  noise  sources  associated  with  the  shunt  resistors. 

The  output  voltage  V  satisfies 

V  =  {♦0/4w)(dA1/dt  +  d<S2/dt).  (22) 

The  quantization  of  the  fluxoid  through  the  outer  superconducting  loop 
yields 


(<5j  -  62)/2lJ  =  (LSK  +  LJ  +  *a)/V 


(23) 


\ 


11 


where  K  is  the  current  through  the  "coupling"  loop.  The  screening  currents 
K  and  J  are  related  via 


l_sd2K/dt2  +  R$dK/dt  +  K/C$  =  R$dJ/dt  +  J/C$  +  dVs/dt 


(24) 


Equations  20-24  are  the  basic  equations  describing  the  double  junction 
SQUID.  For  convenience,  we  re-express  Equations  20-24  in  dimensionless  units 
(see  Table  I).  The  voltage  v  across  the  output  is 


(25) 


where 


=  i  -  j  -  sin  \  +  vnl, 

d62  i  . 

dn  =  2  +  J  '  Sln  (S2  +  vn2 ’ 

=  7i PJ  +  +  2n4>a » 


(26) 

(27) 

(28) 


and 


a  a 

7T^spL  — 2 
s  L  dfr 


6SP 


dk 

do 


+  k 


Bsp  SO 


U  + 


J  +  i, 


(29) 


Note  that  the  noise  current  i  is  related  to  the  dimensionless  Johnson 
noise  voltage  v$  of  the  loss  resistor  R$  vis  i$  =  8sdvs/de. 

Equations  25-29  are  the  basic  equations  describing  the  behavior  of 
the  double  loop  magnetometer.  Note  that  these  equations  reduce  to  the 
conventional  dc  SQUID  equations  in  the  limit  p,  =  0.  In  that  case,  we 
have  seen  that  the  equations  of  motion  can  be  derived  from  a  potential 
u  (Sj^).  However,  in  the  general  case  (t^  f  0),  the  equations  of  motion 
for  the  double  loop  SQUID  cannot  be  derived  from  a  potential  field.  Thus, 
the  thermal  activation  techniques  used  to  solve  the  conventional  dc  SQUID 
cannot  be  applied  to  the  double  loop  structure.  However,  as  we  shall  see 
in  the  next  section,  there  is  a  range  of  parameters  over  which  the 
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TABLE  1 


FLUX  QUANTUM 
APPLIED  FLUX 

JUNCTION  SHUNT  RESISTANCE 
JUNCTION  CRITICAL  CURRENT 
BIAS  CURRENT 

CIRCULATING  CURRENT  (SQUID) 
CIRCULATING  CURRENT  (LOOP) 
OUTPUT  VOLTAGE 
NOISE  VOLTAGE  (SHUNT) 

NOISE  VOLTAGE  (LOSS) 

LOSS  RESISTANCE 
SQUID  SCREENING  FACTOR 
LOOP  SCREENING  FACTOR 
TOTAL  SCREENING  FACTOR 
SQUID  HYSTERESIS  PARAMETER 
SHUNTING  CAPACITANCE 
NOISE  FACTOR 

VOLTAGE  SPECTRAL  DENSITY 
CURRENT  SPECTRAL  DENSITY 

CORRELATION  SPECTRAL  DENSITY 

TRANSFER  FUNCTION 
FLUX  NOISE 
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conventional  and  double  loop  SQUID  characteristics  are  closely  related. 

B.  Double  Loop  Characteristics  Above  Tank  Circuit  Resonance 

First,  the  double  loop  and  conventional  dc  SQUID  equations  are 
identical  for  all  values  of  the  device  parameters,  provided  that  the 
circulating  screening  current  k  =  j  =  constant.  In  that  case,  the 
phase  drops  6j  and  bp  are  constant  and  the  double  loop  SQUID  is  in  the 
zero  voltage  state.  The  screening  currents  j  and  k  are  identical  to 
the  screening  curent  jy  that  would  flow  around  a  conventional  dc  SQUID 
with  total  loop  inductance  Ly  =  L  +  L$  (By  =  6  +  Bg).  Similarly,  the 
dependence  of  the  double  loop  critical  current  on  the  applied  flux  is 
identical  to  the  conventional  SQUID  dependence  for  a  device  with  screen¬ 
ing  factor  By.  The  conventional  dc  SQUID  dependence  of  critical  current 
vs.  applied  flux  is  plotted  in  Figure  10  for  convenience  (from  Ref.  14). 
Note  in  particular  that  the  modulation  depth  Aic  =  i  c  ( s‘>  =  0)  -  i c =  0.5) 
decreases  as  By  increases. 

If  the  double  loop  SQUID  is  biased  above  the  critical  current  in 
the  noise  free  state,  the  phase  differences  and  screening  currents  begin 
to  oscillate  in  time.  For  arbitrary  initial  conditions,  a  transient 
will  be  observed  in  the  screening  currents.  From  Eq.  24,  the  transient 
response  is  given  by  the  resonant  frequency  and  Q  of  the  tank  circuit 
composed  of  L$,  R$  and  Cg  in  series.  The  transient  response  significantly 
complicates  the  numerical  analysis  of  the  device.  A  simplification  occurs, 
however,  if  the  resonant  frequency  of  the  tank  circuit,  =  (lsCs)~2,  lies 
well  below  the  natural  frequency  of  the  junctions,  =  (2tt/4>0)  V  where 
V  is  the  average  voltage  developed  across  the  device  biased  at  some  cur¬ 
rent  i  '■  ic>  In  this  case,  the  Josephson  oscillations  of  the  phase  dif¬ 
ferences  6y  and  62  drive  the  tank  circuit  well  of  resonance.  The 
oscillations  in  j  at  the  Josephson  frequency  and  harmonies  are  essen¬ 
tially  shorted  by  the  capacitance  Cg,  and  the  oscillations  in  the 
screening  current  k  can  be  neglected.  In  terms  of  the  dimensionless 
parameters  B  =  2LI  /<J>  ,  B.  =  2L  I  /<f  ,  B  =  2-ttI  R^C  /*  ,  and  resistance 

U  U  L  U  U  J  O  J  nO 

ratio  p  =  R$/R,  this  corresponds  to  the  condition  <v  •  ""  1.  In 
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addition,  we  restrict  the  discussion  to  a  fairly  high  Q  tank  circuit, 

6$p  <v>  £  1.  This  is  a  physically  reasonable  assumption,  since  a  low 
Q  tank  circuit  would  introduce  excess  Johnson  noise  which  would  tend  to 
degrade  the  device  performance.  The  screening  current  k  can  now  be 
approximated  by  the  time  average,  *k-  =  <j >.  Thus,  Eq.  28  becomes 

-  ^2  =  ni-Sj  +  2n  (4>a  +  By  <j"/2).  (30) 

Equations  25-27  and  Eq.  30  correspond  to  a  set  of  equations  iso¬ 
morphic  to  the  dc  SQUID  equations  with  a  total  applied  flux  increased 
by  the  screening  flux  By  < j  / 2 .  Thus,  the  solution  of  the  conventional 
dc  SQUID  equations  for  screening  factor  8  and  applied  flux  is  an 
approximate  solution  to  the  double  loop  equations  with  screening  factors 
B  and  By  and  applied  flux  <j>a  if  <fa  =  4>y  -  Ry  <j'72.  Note  that  <j>  is  the 
average  circulating  current  for  both  the  conventional  dc  SQUID  solution 
and  for  the  double  loop  solution.  Similarly,  the  average  voltage,  <v>, 
developed  across  the  conventional  dc  SQUID  biased  at  i  and  is  approxi¬ 
mately  equal  to  the  average  voltage  developed  across  the  double  loop  device 
biased  at  i  and  <t>  . 

3 

The  common  solutions  can  be  found  using  a  ioad  line  technique.  We 
illustrate  the  technique  for  a  double  loop  SQUID  with  screening  factor 
B  =  1  and  By  =  50.  We  plot  the  average  circulating  current  <j>  versus 
applied  flux  <j>y  for  a  conventional  dc  SQUID  with  screening  factor  B  =  1 
as  a  function  of  bias  current  in  Figure  H.  Note  that  the  family  of  curves 
are  periodic  in  <t>y  with  period  1  (in  dimensioned  units,  the  periodicity 
is  1’  ).  Only  the  section  0  ij)y  £  1  has  been  shown.  A  load  line,  <j''  = 

(2/B,  )  ( d>T  -  t>  ),  has  also  been  drawn  for  =  0.25,  R.  =  50.  Note  that 

the  load  line  is  not  periodic  in  <f>y.  As  a  result,  a  family  of  parallel 

lines  appears  in  the  collapsed  plot.  Figure  11.  As  the  screening  factor 

By  increases,  the  slope  of  the  load  line  decreases,  and  the  separation 
between  the  parallel  lines  decreases.  As  a  result,  as  By  increases,  the 
load  line  will  intersect  the  curve  <j>  vs.  <)>yfor  fixed  i  bias  at  numerous 
places.  Each  point  of  intersection  represents  a  distinct  approximate 
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solution.  The  corresponding  average  voltage,  <v> ,  can  be  read  off  the 
<v>  vs.  <t>y  characteristics  plotted  in  Figure  12.  In  the  limit,  By  -*  0, 
the  load  line  approaches  the  vertical.  In  this  case,  $>y  =  <j>  and  the 
conventional  dc  SQUID  and  double  loop  solutions  are  identical. 

Noise  free  I-V  characteri sties  for  a  conventional  dc  SQUID  with 
6  =  1  are  plotted  for  various  values  of  applied  flux  in  Figure  4.  All 
the  characteristics  are  single  valued,  and  lie  between  the  characteristics 
labled  4>,  =  0  and  <y  =  0.5.  Noise  free  I-V  characteristics  as  a  function 
of  applied  flux  for  a  double  loop  SQUID  with  3  =  1  and  6,  =  10  are  plotted 
in  Figure  13.  Since  the  average  voltage  at  fixed  bias  current  i  for  the 
double  loop  SQUID  corresponds  to  the  conventional  dc  SQUID  voltage  at  some 
flux  0  <  <j)j  <  0.5,  all  the  double  loop  characteristics  also  lie  within 
the  conventional  SQUID  characteristics  labled  g-y  =  0  and  4>y  =  0.5.  In 
fact,  since  <j>  =  0  at  <>y  =  0  and  4>y  =  0.5,  the  conventional  character¬ 
istics  and  the  double  loop  characteristics  coincide  at  those  values  of 
applied  flux.  However,  as  can  be  seem  from  Figure  13,  a  family  of  hysteretic 
solutions  also  exist.  The  hysteresis  arises  because  the  critical  current 
for  the  double  loop  SQUID  is  identical  to  the  critical  current  for  the 
conventional  SQUID  with  total  screening  factor  By  =  3  +  6|_.  As  By  increases, 
the  critical  current  i  (<j>  )  also  increases  (see  Figure  4b).  However,  as 

C  a 

<v>  increases,  the  double  loop  solutions  are  approximated  by  conventional 
solutions  for  screening  factor  6=1.  These  characteristics  lie  below  the 
characteristics  for  By  =  11.  As  a  result,  the  double  loop  characteristics 
droop.  In  addition,  for  sufficiently  large  values  of  By,  there  are  a  set 
of  metastable  solutions  which  occur  below  the  critical  current.  Another 
set  of  lobed  solutions  occur  above  the  critical  current  (see  Figure  14). 
These  solutions  correspond  to  excess  flux  trapped  in  the  SQUID  loop,  and 
cannot  be  generated  from  the  conventional  solutions  by  smoothly  varying 
the  parameter  By. 

Equations  25-29  were  also  simulated  numerically.  The  equations  were 
integrated  stepwise  in  time  using  a  fourth-order  Adams  Moulton  predictor- 
corrector  method.  The  instantaneous  voltage  v(n)  and  circulating  currents 
j(o)  and  k ( o)  were  computed  at  each  step  interval.  The  initial  conditions 
were  carefully  chosen  so  as  to  minimize  the  transient  response. 
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The  average  voltage  was  determined  by  computing  the  period  of  oscillation 
of  the  phases  and  ^  I-V  characteristics  determined  by  this  method 
are  in  excellent  agreement  with  the  approximate  hysteretic  solutions. 
Appropriate  initial  conditions  could  not  be  determined  for  the  negative 
resistance  region  of  the  hysteretic  characteristics,  nor  for  the  meta¬ 
stable  lobes. 

C.  Noise  Free  Characteristics  for  Arbitrary  Parameter  Values 

Noise  free  characteristics  for  arbitrary  values  of  the  circuit  para¬ 
meters  are  generated  by  direct  numerical  simulation  of  the  basic  equations. 

In  the  general  case,  no  simple  relationship  exists  between  the  double 
loop  characteristics  and  the  conventional  dc  SQUID  characteristics.  We 
are  interested  in  calculating  the  characteristics  for  values  of  the  para¬ 
meters  which  enhance  the  coupling  properties  of  the  device  without  degrading 
the  signal-to-noise  ratio.  A  reasonable  quess  would  be  to  choose  values 
for  the  resonant  coupling  circuit  consisting  of  the  coupling  loop  Ls, 
shorting  capacitance  C$  and  loss  Rs  such  that  the  resonant  frequency  of 
the  coupling  circit  is  much  less  than  the  Josephson  frequency  of  the  junc¬ 
tions,  <<  n«j.  As  demonstrated  in  Section  III-B,  the  I-V  characteristics 
in  that  case  become  multivalued  and  show  considerable  hysteresis.  Unfor¬ 
tunately,  hysteresis  in  the  I-V  characteristics  of  a  conventional  dc  SQUID 
near  the  zero  voltage  state  degrades  the  performance  of  the  device.  To 
avoid  this  problem  in  the  double  loop  SQUID,  we  must  relax  the  constraint 

(,'t  <<  V 

In  this  section,  a  family  of  I-V  characteristics  is  computed  for 
double  loop  devices  with  variable  resonant  frequencies  If  the  coupling 
loop  inductance  is  fixed,  the  resonant  frequency  <>'t  is  a  function  of 
the  dimensionless  shorting  capacitance  B$,  uj^  =  ( n 0^ ) ~ ^ •  We  take  the 
resistance  ratio  p  =  R$/R  <<  1  to  reduce  the  Johnson  noise  generated  in 
the  loss  resistance  R$  (the  noise  free  characteristics  are  essentially 
independent  of  p  in  this  limit). 

The  I-V  characteristics  generated  by  nimerical  simulation  are  plotted 
in  Fig.  15  as  a  function  of  i,^.  The  screening  factors  are  i-;  =  1.0,  =  10.0 
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and  the  applied  flux  <f>  =  0.25.  For  comparison,  the  dimensionless  Josphson 

a 

frequency  <oj  =  <v  •.  Thus,  the  curve  Tabled  3$  _  10,  c..  0.06  corresponds 

to  the  hysteretic  approximate  solution  for  ^  <<  Wj  found  in  Section  III-B. 

Note  that  the  negative  resistance  region  and  metastable  lobes  were  not 
generated  in  the  numerical  simulation.  As  nk  is  increased,  the  hysteresis 
near  the  zero  voltage  state  increases.  In  fact,  the  characteristic  plotted 
in  Fig.  13  for  .,k  =  0.15  lies  below  the  conventional  dc  SQUID  characteristic 
for  3  =  1,  .-v  =  0.5.  In  this  case,  the  double  loop  behavior  clearly  cannot 
be  approximated  by  the  conventional  dc  SQUID  behavior.  However,  at  bias 
currents  such  that  u>,  =  <v  ■  N>  0.15,  the  characteristic  approaches  the 
approximate  solution  found  in  Section  III-B. 

The  I-V  characteristic  for  =  0.3  shows  a  qualitatively  different 
behavior.  In  this  case,  the  tank  circuit  resonant  frequency,  uv,  corresponds 
to  the  average  voltage  near  the  bottom  of  the  droop  in  the  approximate 
solution.  As  a  result,  for  :v  £  w^.,  the  Josephson  oscillations  are  not 
effectively  shorted  out  by  the  capacitance  C$.  The  device  behaves  like 
a  conventional  dc  SQUID  with  3  =  11,  <}>  =  0.25.  No  hysteresis  is  observed 

^  N, 

near  the  zero  voltage  state.  In  the  region  <v  ^  u>t,  the  double  loop 
characteristic  begins  to  deviate  from  the  3  =  11  dc  SQUID  characteristic. 
Furthermore,  the  numerical  simulation  becomes  somewhat  unstable. 

The  character  of  this  instability  can  be  seen  in  Fig.  16.  The  instantaenous 
voltage  v  and  screening  currents  k  and  j  are  plotted  as  a  function  of  time, 

0,  for  a  double  loop  SQUID  with  3  =1,  6L  =  10  and  3  =  0.3537.  At  bias  cur¬ 
rent  i  =2.0,  (uij  <  )  the  Josephson  frequency  u>j  <  i^,  and  the  voltage  and 

screening  currents  are  periodic  functions  of  time  as  expected.  However,  if 
the  bias  current  is  increased  to  i  =  2.25,  the  Josephson  frequency  approaches 
the  resonant  frequency  of  the  coupling  loop.  At  this  point,  the  voltage  and 
circulating  currents  are  no  longer  periodic  functions  of  time.  If  the  bias 
current  is  then  increased  to  i  =  2.4,  utj  N  ■  and  the  periodic  relationships 
between  v,  j  and  k  and  the  time  o  return.  This  kind  of  behavior  is  character¬ 
istic  of  chaotic  systems.  In  this  case,  the  equations  of  motion  can  be  reduced 
to  a  set  of  coupled  third  order  non-linear  differential  equations  describing 


a  system  of  two  lightly  damped  non-linear  oscillators  with  nearly  equal 

characteristic  frequencies.  Chaotic  behavior  in  systems  of  this  nature  has 
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been  described  in  detail.  ’  From  a  practical  standpoint,  we  expect  chaotic 
behavior  to  be  accompanied  by  a  continuous  low  frequency  voltage  power  spectral 
density,  even  in  the  absence  of  Johnson  noise  sources.  Although  chaotic 
behavior  in  Josephson  junction  devices  is  a  subject  of  considerable  interest, 
a  careful  investigation  of  this  behavior  in  the  double  loop  device  is  outside 
the  scope  of  the  present  effort. 

At  bias  currents  i  c  2.2,  instabilities  in  the  calculation  of  the 
I-V  characteri stic  labled  .i.  =  0.3  are  no  longer  observed.  In  this  region, 
a  hysteresis  loop  is  observed.  Note  that  both  branches  of  the  loop  are 
numerically  stable.  At  bias  currents  i  >  3.0,  the  characteristic  approaches 
the  approximate  solution  for  the  limit  ov  «  wj.  In  this  region,  both 
solutions  are  nearly  equal  to  the  conventional  dc  SQUID  solution  for  *  =  0.5, 

a 

3  =  1. 

The  final  characteristic  plotted  in  Fig.  15  corresponds  to  the  limit 
11 1  N''  ‘J’  ^or  f re(luenc''es  li;j  within  the  hysteretic  region  in  the  approximate 
solutions  (i.e.,  _  0.6).  In  this  case,  the  I-V  characteristic  is  nearly 

identical  to  the  conventional  dc  SQUID  characteristic  for  p  =  11,  <t  =  0.25. 

a 

Furthermore,  no  instabilities  in  the  simulation  were  observed.  This  is 
probably  because  the  non-linearities  associated  with  the  Josephson  junctions 
are  small  in  the  region  o>t  n,  mj.  As  a  result,  chaotic  effects  are  less 
likely. 

In  conclusion,  the  only  characteristics  in  Fig.  15  which  are  non- 
hysteretic  satisfy  o,  0.6.  In  this  case,  the  shorting  capacitance  is 
essentially  irrelevant,  and  the  characteristics  reduce  to  the  conventional  dc 
SQUID  characteristic  with  1^.  =  P  +  fi^.  As  is  decreased,  the  I-V  character¬ 
istics  begin  to  deviate  from  the  conventional  dc  SQUID  characteristics. 

In  particular,  a  hysteresis  loop  is  observed  at  bias  currents  i  s  ic> 

For  sufficiently  low  values  of  u^,  the  hysteresis  occurs  near  the  zero 
voltage  state. 


D.  Numerical  Techniques  for  Double  Loop  Characteristics  with  Noise 

In  this  section,  the  numerical  techniques  used  to  compute  device  char¬ 
acteristics  in  the  presence  of  noise  are  discussed.  The  dominant  noise 
sources  within  the  device  are  assumed  to  be  the  Johnson  noise  generated 
by  the  junction  shunt  resistances  and  the  tank  circuit  resistance,  R  .  The 
effect  of  low  frequency  noise  sources  will  not  be  discussed  in  this  section. 

We  assume  that  the  device  parameters  and  operating  point  are  chosen  such 
that  the  power  spectral  density  of  the  voltage  noise  source  associated  with 
each  resistance  is  white  at  all  relevant  frequencies.  For  a  device  generated 
at  temperature  T,  the  corresponding  average  voltage  across  the  device  must 
satisfy  V  =  <v>I0R  «  27rkBT/s'o.  The  voltage  noise  sources  are  uncorrelated, 
with  shunt  power  spectral  density  =  4kgTRs  and  loss  power  spectral  density 
S$  =  4kgTRs-  In  dimensionless  units,  Sn  =  4r  and  Ss  =  41’n. 

In  the  conventional  dc  SQUID,  the  equations  of  motion  written  in  dimen¬ 
sionless  form  are  independent  of  the  shunt  resistance,  R.  As  a  result, 
the  energy  factor  scales  with  the  shunt  resistance,  o  1/R.  The  value 
of  R  is  limited  by  the  constraint  £c  =  2nl  R^C  <  1.  The  dimensionless  equa¬ 
tions  of  motion  for  the  double  loop  SQUID  are  also  independent  of  R.  How¬ 
ever,  the  equations  are  a  function  of  the  resistance  ratio  r  =  Rs/R-  As 
we  have  seen  in  Section  III-C,  the  noise  free  I-V  characteristics  are  rela¬ 
tively  independent  of  t>,  provided  P$p<v''  <  1.  In  the  noise  rounded  v> 
however,  any  noise  generated  by  the  loss  resistor  will  tend  to  degrade  the 
resolution  of  the  device.  Thus,  we  are  interested  in  computing  the  character¬ 
istics  for  device  parameters  such  that  the  noise  currents  generated  by  the 
loss  resistance  are  small  compared  to  the  noise  currents  generated  by  the 
junction  shunt  resistance. 

An  indication  of  the  relative  importance  of  the  noise  currents  generated 
in  the  device  is  obtained  by  taking  the  limit  I0  =  0.  In  this  limit,  the 
shunted  junctions  are  replaced  by  the  shunt  resistances  (see  Figure  17). 

The  device  now  consists  of  ordinary  linear  circuit  elements  with  uncor¬ 
related  white  voltage  noise  sources  in  series  with  the  resistances  R  and 
Rs<  The  voltage  noise  source  V$  in  series  with  R$  generated  a  current  Is 
through  the  shorting  capacitance  Cs.  This  noise  source  affects  the  device 
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in  two  ways.  First,  noise  components  near  the  Josephson  frequency  <.>j  will 
contribute  to  the  noise  rounding  of  the  I-V  characteristics.  In  addition, 
the  noise  will  be  mixed  down  by  the  non-linear  interaction  in  the  junctions 
and  appear  as  a  contribution  to  the  low  frequency  voltage  and  current  noise. 

In  addition,  the  component  of  I$  near  the  tank  circuit  resonance,  .t ,  can 
link  noise  flux  into  the  device  through  the  inductances  L  and  Ls<  Since 
the  I-V  characteristics  and  device  noise  spectral  densities  are  a  function 
of  the  quasistatic  applied  flux  ,  the  device  characteristics  at  frequencies 
m  <<  m.  can  be  approximated  by  a  weighted  average  of  the  characteristics 
over  the  variance  in  the  flux  produced  by  the  noise  in  the  loss  resistor 
R  .  If  the  flux  resolution  is  a  sharp  function  of  the  applied  flux,  the 
average  resolution  may  be  reduced  considerably. 

As  an  example,  consider  a  double  loop  SQUID  with  ft  =  1,  i- ^ =  L$/L  1. 
Let  the  resonance  frequency  of  the  tank  circuit  aij.  <<  m,0  =  The 

noise  current  I$  =  vs/^j  where 

Zy  =  R$  [1  +  2fv>QZ(.,/...t)Z]  +  jRs  CvQ(4  +  AlT/R2)(./.t)  -  Q(..t/J]  (31) 


and 

y  =  (4  +  m2L2/R2)  1  (32) 

Case  A:  i.^Ly  s  a>tLs  vv-  2R.  In  dimensionless  units,  pQ  2.  This 
condition  is  equivalent  to  3|_m^/i.y0  1.  At  frequencies  ^  y  ~ 

R2/iA2,  and  Zy  R$  (1  +  2/p)  +  jR$Q/PL.  The  fraction  of  the  current  noise 
from  R^  flowing  around  the  SQUID  loop  is  Sjs  -  4kgTR$/ j Zy | 2  '  4kgT/Rs 
[(1  +  2/n)2  +  Q2/p2].  By  comparison,  the  current  noise  from  the  shunt  re¬ 
sistances  flowing  in  the  SQUID  loop  is  Sj^  4kgT(2R)/R2[(2  +  p)2  +  ( TT •  /'  j0  " 
Qp(„t/.,)2].  Thus  the  noise  from  the  loss  resistor  can  be  neglected  if 
(2/o)[(2  +  p)2]  v  (2  +  p)2  +  [ n (u)/u.j0 )  -  Qp  (u't/i')]2  (note  pQ/i-y  «  1). 

Now  Qp(a>t/,v)  ~  Qp  (mt/wJo)  2ftL(.,»t/u!Jo)2.  Thus,  we  need  ftL  (wt/uiJo)Z 
p  «  4,  p  <<  1.  If  p  is  small  enough  to  satisfy  both  of  these  conditions, 
the  effects  of  the  noise  from  R$  on  the  junctions  can  be  neglected.  At 


21 


•'***'*- 


frequencies  near  the  resonant  frequency. 
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For 


:<  1,  S 


0  *  4kBT/RsQ^. 

^  _ _  _  i _ i 


The  noise 


4kgTL  /RCQ  over  a  bandwidth  Af  =  i4/2t!Q. 


ZT  ~  Rs  (1  +  2/p)  +  jRsQ. 

flux  spectral  density  is  S, 

2  2  u  J  q 

The  variance  .'4  /$Q  '  {r/2n )  (1/oQ  )  <<  1.  Thus  the  noise  flux 

through  the  SQUID  loop  can  be  neglected.  The  noise  current  through  the 
coupling  loop  S ^  Sj5  (2R/:.i(.L(.  produces  a  noise  flux  of  S,c  = 

LJSKs  (Ls/L)2(2RA,tLs)2  S.s.  The  variance  A42/.^  (r/2n)  (^Q/  ,t ) 

( 1  / oQ 3 )  =  (i716TT)(„),  /at)4(o2/a3)-  Since  J,  '>  1  and  r,  (•*/■•!„)  1, 


Jo'  t'  '  '  L  '  v't  p|_  1  a"u  lLv‘t/,‘Jo' 

this  condition  also  places  an  upper  limit  on  p.  Note  that  the  amount 

of  variance  in  flux  at  ai*  which  can  be  tolerated  without  serious  loss 
of  resolution  depends  on  the  detailed  dependence  of  the  energy  resolu¬ 
tion  on  applied  flux. 

Case  B:  a^Ly  '  Myl_s  ■-<  2R.  In  dimensionless  units  pQ  <<  2  and 
bL (‘t /  ‘Jq )  <<-  1*  Near  the  resonance  frequency,  \  0.25  and  Zy  ? 

(1  +  /2o)  +  j RSQ 3K)3/4i 1-^.  The  flux  noise  generated  through  the 

coupling  loop  is  L2  Svs/|ZT|2  4kBTRsL2/R2(I+,  Q2)2.  Since 

3 


1,  the  averaging  over  the  characteristics  produced  by  V  can  be 

,2 


!  "  »  -  -  •  -  •  -  j  •  j  —  ■  —  •  -  - - -  r  - - - —  j  ’  c 

L  2  b 

large  unless  pQ  1.  Since  pQ  «  2,  this  places  a  constraint  on  both 
p  and  Q,  1/Q2  «  p  <-•;  2/Q  and  Q  ■  1. 

In  conclusion,  the  effect  of  the  noise  source  Vs  on  the  device 
characteristics  can  be  ignored  for  a  wide  range  of  device  parameters. 
For  a  double  loop  SQUID  with  resonant  frequency  ^  <  u'jQ  and  coupling 


loop  inductance  L 


L,  the  requirement  is  typically  R$  <<  R. 


We  assume  that  the  double  loop  SQUID  parameters  are  chosen  such 
that  the  noise  in  V  can  be  ignored.  In  this  case,  the  device  character¬ 
istics  are  computed  as  follows.  Equations  25-29  are  integrated 
stepwise  in  time  on  a  computer.  The  noise  sources  in 

the  shunt  resistors,  v  £  and  v  £,  are  approximated  as  follows.  The  Johnson 
noise  terms  are  considered  to  consist  of  a  chain  of  pulses  of  equal  dura¬ 
tion  ad  and  random  amplitude  v^.  The  amplitudes  {v^l  should  be  chosen 
such  that  the  voltaqe  noise  power  spectral  density  is  white,  with  value 
Sn(f)  =  4r  for  I  =  ^TTkgT/ IQ  :-0  (in  dimensioned  units,  SN(F)  =  4 k^TR ) . 
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To  facilitate  the  computation,  the  noise  source  is  artificially  smoothed 
by  taking  the  Courier  transform  of  the  voltage  amplitudes  =  a  exp  [ i , 
where’  is  a  pseudorandom  number  uniformly  distributed  in  the  range 
(0,  2v).  Thus  the  spectral  density  S  (f)  =  2 T  at  all  values  of  f  for 
each  individual  run.  The  total  nunber  of  pulses  in  a  run,  N,  and  the 
pulse  duration,  AO,  are  chosen  such  that  the  minimum  frequency  component 
fm-n  =  1/NA0  <<  where  m^is  the  tank  circuit  resonant  frequency.  In  addi¬ 
tion,  the  maximum  frequency  component  must  satisfy  f  -  1/2' A  v/here 

.  is  the  fundamental  frequency  of  the  Josephson  oscillations.  Since  we 

ary;  interested  in  computing  the  characteristics  for  the  case  these 

t  J 

conditions  imply  that  K  must  be  a  large  number.  Typical  values  are 

fmin  :‘t/20’  t  *  V10*  and  N  V  ?048- 

The  differential  equations,  Eqs.  11-15,  are  then  integrated  step¬ 
wise  in  time  using  a  fourth-order,  Adams  Moulton  Predictor-corrector 
method.  The  instantaneous  voltage  v(0)  and  circulating  currents ,  j(s) 
and  k  ( 0 )  are  computed  at  each  step  interval.  Noise  rounded  i-v  character¬ 
istics  are  calculated  by  computing  the  averaqe  voltage  <v>  over  typically 
AO  sets  of  2048  intervals.  The  initial  values  for  the  integration  are 
carefully  chosen  so  that  no  transient  response  is  apparent  in  the  average 
voltage.  Noise  power  spectral  densities  Sy ,  S-,  $k,  Syj  and  Svk  are 
determined  from  the  average  of  the  products  of  t.ho  Fast  Fourier 
Transform  of  the  appropriate  variables  for  a  single  run.  We  estimate 
the  average  voltage  to  be  accurate  to  about  1*  and  the  noise  power  spectral 
densities  to  be  accurate  to  about  15X 


E.  NOISE  ROUNDED  CHARACTERISTICS  FOR  THE  DOUBLE  LOOP  SQUID 

In  this  section  we  calculate  noise  rounded  characteristics  for  two 

particular  devices.  First,  we  consider  a  double  loop  SQUID  with  P  =  1, 

=  1  and  n.Q  «  Mjo.  A  device  which  is  similar  to  this  structure  is 

presently  under  construction  by  Ho  Jung  Paik  and  co-workers.  The  device 

has  been  fabricated  in  a  toroidal  geometry  with  point  contacts  as  the 
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Josephson  junctions.  In  this  case,  the  coupling  properties  ate  enhanced 
by  the  properties  of  the  toroidal  geometry,  rather  than  by  an  increase 
in  the  coupling  loop  6^.  Next,  we  consider  a  double  loop  SQUID  with 
8  =  1 ,  8.  =  10.  In  this  device  the  coupling  loop  inductance  has  been 
increased  by  an  order  of  magnitude.  Significant  deviation  from  the  con¬ 
ventional  dc  SQUID  behavior  are  expected  in  this  case. 


1.  Double  Loop  SQUID  with  8  =  1 ,  (%|_  =  I 

The  device  parameters  are  chosen  such  that  the  resonant  frequency 

2 

u)  «  u>,„,  and  the  resistance  ratio  p  < <•  1.  We  take  l-v  =  2-1  R  C/tp  =  100 

O  00  S  o  s  u 

and  p  =  R  /R  =  0.01.  Calculations  are  performed  for  a  moderately  noise 
rounded  device,  F  =  0.01.  The  resonant  circuit  Q  =  (1/p  )(*8l/Ps)5  =  18, 
and  the  ratio  0^/wjo)  =  0.056.  Note  that  oQ  ~  0.2  and  the  device  is 
near  the  case  B  limit.  In  particular,  1/Q^  ?  0.003  <  e  <.  2/Q  ?  0.1.  The 
ratio  )  and  the  resonant  circuit  Q  are  chosen  to  facilitate  the 

nunerical  calculation.  Note  that  the  calculation  has  been  performed  for 
a  relatively  low  value  of  Q.  In  general,  a  transient  will  be  generated 
in  the  circuit  at  the  beginning  of  the  numerical  simulation.  The  size 
of  the  transient  is  determined  by  the  initial  values  of  the  phases  and 
circulating  currents.  The  transient  decay  time  is  determined  by  the 
Q  of  the  tank  circuit.  In  the  case  .vt  <<  ,  Q  1,  the  integration 

must  be  carried  out  for  a  large  number  of  Josephson  periods  before  useful 
data  is  generated.  To  avoid  this  problem,  we  have  performed  the  calcula¬ 
tion  for  a  relatively  low  value  of  Q.  We  expect  the  results  to  be  rela¬ 
tively  independent  of  Q,  provided  Q  1,  p  <•;  1  and  <  .)jQ  <<  1. 
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A  family  of  characteristics  for  a  double  loop  SQUID  with  the  para¬ 
meters  listed  above  is  plotted  in  Figure  IB.  For  comparison,  we  also 
plot  noise  rounded  characteristics  for  a  conventional  dc  SQUID  with  the 
same  total  loop  inductance  By  =  B  +  py  =  2.  In  principle,  the  input 
signal  can  be  coupled  into  the  entire  superconducting  loop.  Thus  the 
input  inductance  of  the  device  is  given  by  By,  rather  than  by  By.  As 
a  result,  the  energy  resolution  of  the  double  loop  SQUID  should  be  com¬ 
pared  to  the  resolution  of  a  conventional  dc  SQUID  with  screening  factor 

i'T- 

The  double  loop  I-V  characteristics  are  plotted  in  Figure  18a.  The 
noise  free  characteristic  labled  r  =  0  is  hysteretic,  with  critical  cur¬ 
rent  i  =  1.75.  Solutions  in  the  negative  resistance  region  of  the  character¬ 
istic  could  not  be  generated  numerically.  Note  that  the  noise  rounded 
characteristic  labled  1'  =  0.01  lies  below  the  noise  free  curve.  In  the 
conventional  dc  SQUID,  the  opposite  behavior  is  observed.  Also  note 
that  the  double  loop  I-V  characteristic  is  not  well  approximated  by  the 
noise  rounded  dc  SQUID  characteristic  for  ,y  =2. 

The  forward  transfer  function,  >  <v  /  •  j>  ,  is  plotted  as  a  function 

d 

of  bias  current  for  4>a  =  0.275  in  Figure  18b.  At  bias  currents  i  A  2.0, 
the  double  loop  transfer  function  for  f  =  0.01  is  identical  to  the  noise 
free  value.  The  noise  free  transfer  function  becomes  infinite  as  the 
characteristic  approaches  the  negative  resistance  region  (i  o.  1.7).  The 
noise  rounded  characteristic  is  well  behaved  in  this  region.  The  noise 
rounded  double  loop  transfer  function  is  not  well  approximated  by  the 
conventional  noise  rounded  SQUID  characteristic  for  By  =  2.  The  low 
frequency  voltage  and  current  noise  spectral  densities  are  plotted  in 
Figure  18(c)-(g).  Curves  are  plotted  for  a  double  loop  SQUID  with 
<!>  =  0.25  and  <p  =  0.3,  and  for  the  conventional  dc  SQUID  at  a  =  0.25. 

o  a  a 

The  low  frequency  voltage  noise  spectral  density  S°  peaks  up  near  the 
critical  current  i  1.75.  At  this  bias  current,  the  double  loop  value 
is  more  than  twice  the  conventional  dc  SQUID  value  for  By  =  2.  At  bias 
currents  i  ic,  the  voltage  noise  spectral  densities  approach  the 
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spectral  density  of  the  parallel  shunt  resistances,  S°  -  2|'.  The  low 

frequency  current  noise  spectral  densities  S°  and  S°  plotted  in  Fig. 

J  K 

18  (d)-(e)  also  peak  near  i  .  However,  at  this  bias  current,  the  con¬ 
ventional  Pj  =  2  SQUID  and  double  loop  spectral  densities  are  nearly 
equal.  Since  the  conventional  SQUID  current  noise  tends  to  decrease 
for  increasing  screening  factor  (see  Appendix  I),  the  double  loop  current 
noise  is  actually  depressed  below  the  conventional  value  for  SQUID  ,1-  =  1. 

At  bias  currents  i  '•>  ic»  the  current  noise  spectral  densities  approach 

the  value  for  the  series  combination  of  the  shunt  resistances  S° ,  S°  -  2;'. 

J  K 

The  voltage  and  current  noise  spectral  densities  are  strongly  correlated 

near  i  (Fig.  15  (f)-(g)).  At  bias  currents  i  >-■  i  ,  the  noise  sources 
^  c 

are  essentially  uncorrelated. 

The  dimensionless  low  frequency  flux  noise  r°.  =  (S°/2r)/(3<v  /3t  ) \ 

v  v  a 

is  plotted  as  a  function  of  bias  current  in  Fig.  18(h).  From  Table  1,  the 

dimensioned  low  frequency  flux  noise  S°  =  <'°  (2kDTR)/(I  R/tf  )2  = 
o  2  2  '•’•(Boo 

r,t)^kBTi/IoR^o'  This  is  the  ecluivalent  flux  noise  at  the  input  which 

would  produce  the  voltage  noise  observed  at  the  device  output.  Note 
that  the  minimum  dimensionless  flux  noise  for  the  double  loop  SQUID  is 
almost  identical  to  the  minimum  for  the  conventional  dc  SQUID  with  the 
same  total  loop  inductance.  Furthermore,  since  the  dimensioned  flux 
noise  is  a  function  only  of  the  junction  parameters  I  ,  R  and  T,  and 
the  dimensionless  factor  f,.*  the  dimensioned  flux  noise  is  also  the  same 
in  both  cases.  Thus,  although  the  two  devices  has  dissimilar  I-V  character¬ 
istics,  forward  transfer  functions  and  voltage  spectral  densities,  the 
effective  input  flux  noises  are  the  same.  Also  note  that  in  the  conventional 
SQUID,  the  minimum  flux  resolution,  the  maximum  voltage  noise  spectral 
density  and  the  maximum  forward  transfer  function  all  occur  at  approximately 
the  same  values  of  bias  current,  I  %  1.75  I  .  For  the  double  loop  SQUID, 
however,  the  minimum  flux  noise  occurs  at  bias  current  i  >■  1. 9-2.0.  This 
is  the  value  at  which  the  noise  rounded  average  voltage  began  to  decrease 
below  the  noise  free  values  (i.e.,  the  hysteresis  in  the  I-V  curve  began 
to  affect  the  device  performance).  Thus,  in  operatina  this  device,  it 
is  important  to  choose  that  value  of  bias  current  which  optimizes  the 


flux  resolution  (i.e.,  I  =  1.9  I  )  and  not  the  value  which  optimizes 
the  forward  transfer  function  (i.e.,  I  =  1.7  IQ). 

The  magnetic  field  resolution  of  the  double  loop  SQUID  is  proportional 
to  the  energy  factor  S£.  If  the  input  signal  is  coupled  to  the  entire 
superconducting  loop  inductance  L_,  the  energy  factor  is  =  S°/2Ly. 

As  we  have  seen,  the  flux  noise  S°  for  a  double  loop  SQUID  with  :•  =  1 , 

R.  =  1  is  nearly  identical  to  that  of  a  conventional  dc  SQUID  with  the 
same  total  loop  inductance.  Thus  the  energy  factor  and  input  inductance 
of  the  two  devices  are  almost  the  same.  As  a  result,  the  magnetic  field 
resolution  of  the  two  devices  used  directly  as  sensors  will  be  similar. 

The  main  different  is  that  the  double  loop  SQUID  structure  can  be  easily 
realized  in  a  toroidal  geometry.  Thus  if  the  device  is  coupled  through 
a  mutual  inductance  to  a  coupling  coil  in  series  with  a  pick-up  loop, 
the  superior  coupling  properties  of  the  toroidal  geometry  may  be  exploited 
to  improve  the  magnetic  field  resolution. 
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2.  Double  Loop  SQUID  with  (  =1,  = 

Results  for  the  double  loop  SQUID  with  SQUID  loop  :-=  1  and  coupling 
loop  3^  =  10  are  plotted  in  Fig.  19  (a)-(h).  Calculations  are  performed 
for  p  =  0.1  and  P  =  10.  For  these  values,  the  resonant  circuit  Q  =18 
and  the  ratio  =  0.056  as  in  the  previous  case.  The  total  super¬ 

conducting  loop  inductance  Py  =  11.  The  critical  current  in  the  noise  free 
case  is  i  =  1.93  at  .|>  =  0.25.  As  a  result,  the  noise  free  double  loop 

C  cl 

characteristic  plotted  in  Fig.  19  (a)  shows  substantial  hysteresis.  Cal¬ 
culations  are  performed  for  a  relatively  low  value  of  the  noise  figure, 

F  =  0.001,  to  illustrate  the  effects  of  this  hysteresis  on  the  noise  rounded 
characteristics . 

The  noise  rounded  I-V  characteristic  is  plotted  in  Fig.  19(a).  The 
characteristic  is  almost  identical  to  the  noise  free  double  loop  characteristic 
for  i  >  1.9.  Below  this  value,  the  noise  rounded  characteristic  switches 


rapidly  to  the  zero  voltage  state.  At  this  value  of  bias  current,  the  double- 
loop  noise  free  characteristic  was  generated  from  the  conventional  dc  SQUID 
characteristic  for  3  =  1,  and  effective  flux  ~  0.4.  The  critical  current 
of  a  conventional  SQUID  with  B=  1,  .J>  =  0.4,  is  i  =  1.26.  Thus  the  double 

a  C 

loop  noise  free  behavior  in  the  region  i  ^  1.9  is  similar  to  a  conventional 
SQUID  biased  well  above  its  critical  current.  The  noise  rounding  of  the 
I-V  characteristics  of  a  conventional  dc  SQUID  is  negligible  for  F  =  0.001 


in  this  region. 

The  similarity  between  the  double  loop  characteristics  for  }  =  0.25, 

i  ^  1.9  and  the  conventional  dc  SQUID  characteristics  for  P  =  1,  ;  =  0.4 
is  also  apparent  in  the  plots  of  the  double  loop  forward  transfer  function, 
’)<v  /;vi'  ,  and  voltage  noise  spectral  density  S°/2r  (Fig.  19  (b)-(c)).  The 
transfer  function  for  r  =  0.001  is  nearly  identical  to  the  noise  free  value 
for  i  •  1.9  as  expected  for  a  device  biased  well  above  its  critical  current. 
The  voltage  spectral  density  also  rapidly  approaches  the  value  for  a  conven¬ 
tional  device  biased  well  above  the  critical  current,  S°  -  2i\  at  bias  cur¬ 
rents  i  1.9.  For  comparison,  the  voltage  noise  spectral  density  of  a 
conventional  dc  SQUID  with  P,  =  11  is  S°/2r  -  70.  This  discrepancy  between 


the  behavior  of  the  double  loop  SQUID  with  Py  =  11  and  the  conventional  SQUID 
with  the  same  loop  inductance,  ft  =  11,  is  even  more  striking  than  in  Section 
III.E.l.  Below  i  =  1.9,  the  voltage  noise  increased  dramatically  as  the 
SQUID  switches  to  the  zero  voltage  state,  then  falls  to  S°/2r  <  10’*\ 

The  low  frequency  double  loop  SQUID  circulating  current  noise  spectral 
densities,  S^/2r  and  S?/2;',  are  similar  (Fig. 19  (d)-(e)).  Both  values  are 

A 

depressed  below  the  value  for  the  isolated  shunt  resistances,  S^  =  2r,  at 
1.9  <  i  £  2.1.  As  the  SQUID  switches  to  the  zero  voltage  state,  the  current 
noises  increase  slightly.  However,  the  increase  is  significantly  less  than 
that  observed  for  the  voltage  noise,  because  the  device  is  switching  between 
states  of  similar  circulating  current.  The  low  frequency  noise  in  the  current 
j  and  k  are  moderately  correlated  with  the  voltage  noise  (Fig.  19  (f)-(g)). 

The  double  loop  SQUID  flux  noise  spectral  density,  =  (S°/2r)/(9-v  /Vg) 
is  plotted  as  a  function  of  bias  current  for  a  =  0.275  in  Fig.  19  (h).  Above 
i  =  1.9,  the  flux  noise  is  well  approximated  by  the  conventional  dc  SQUID 
value  for  a  SQUID  with  P  =  1,  =  0.4.  The  flux  noise  for  this  device  reaches 

a 

a  minimum  near  i  =  1.26.  However,  the  double  loop  flux  noise  increases 
rapidly  at  i  <  1.9  due  to  the  hysteresis  in  the  I-V  characteristic.  As  a 
result,  the  minimun  double  loop  flux  resolution,  f™in-  90,  is  similar  to  the 

vj) 

minimum  value  for  a  conventional  dc  SQUID  with  P  =  11,  :  =  0.25,  v10'  100. 

Thus,  although  the  device  characteristics  are  similar  to  a  conventional  SQUID 
with  ft  =  1 ,  <j>  =  0.4  for  i  1.9,  the  minimum  flux  resolution  is  limited 

a 

to  a  value  determined  by  the  total  loop  inductance  fy  =  11. 


IV.  ANALYSIS  OF  CONVENTIONAL  DC  SQUID  MAGNETOMETER  SYSTEMS 

In  this  section  a  model  is  developed  for  conventional  dc  SQUID  magnet¬ 
ometer  systems.  Two  examples  are  given.  First,  the  dc  SQUID  is  coupled 
inductively  to  a  superconducting  input  circuit  containing  a  pick-up  coil 
and  an  input  coil  in  series.  Next,  the  effect  of  including  a  resistance  in 
series  with  the  input  inductance  is  investigated.  This  calculation  is  a 
aeneral ization  of  the  model  for  the  dc  SQUID  magnetometer  described  by 
Clarke,  Tesche  and  Giffard  (CTG)?^ 

A.  DC  SQUID  with  Superconducting  Input  Circuit 

A  conventional  dc  SQUID  coupled  inductively  to  a  superconducting  input 

circuit  is  shown  in  Fig.  20  (a).  The  external  quasistatic  flux  *  is  applied 

6 

to  a  pick-up  inductance  L  in  series  with  an  input  inductance  L-.  The  mutual 

P  2  2 

inductance  between  the  input  inductance  and  the  SQUID  loop  is  M  =  n  LL^ . 

The  effects  of  stray  capacitance  within  the  input  circuit  and  between  the 

input  circuit  and  the  SQUID  loop  are  not  included  in  this  model. 

Flux  is  quantized  in  the  input  circuit.  As  a  result,  the  circulating 

current,  J,  in  the  SQUID  loop  generates  a  superconducting  screening  current, 

H,  in  the  input  circuit  satisfying 

(Li  +  Lp)  H  +  MJ  +  <J'a  =  0.  (33) 

Quantization  of  the  fluxoid  in  the  SQUID  loop  implies 

D  =  LJ  +  MH  (34) 

where  D  =  ($o/2n)  (^  -  ?)■  Combining  Eqn.  33  and  34, 

D  =  LeJ  +  4>e,  (35) 

where 

Le  «  L ( 1  -  nr)  (36) 
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(37) 


!e  =  .xr(L/Li)* 


and 


r  =  Li/(Li  +  Lp).  (38) 

Eqn.  35  is  the  analog  of  the  fluxoid  quantization  condition  for  an  isolated 
dc  SQUID.  The  equations  of  motion  for  the  phase  drops  and  <S^  are  the  usual 
junction  relations,  Eqn.  26  and  27.  The  voltage  across  the  SQUID  output  satisfies 
Eqn.  25.  As  a  result,  the  equations  for  the  superconducting  dc  SQUID  magneto¬ 
meter  are  identical  to  those  of  an  isolated  dc  SQUID  with  loop  inductance  L 

1 4  e 

and  applied  flux  7 

The  superconducting  dc  SQUID  magnetometer  may  be  modeled  as  a  conventional 
dc  SQUID  with  input  inductance  Lp.  The  device  forward  transfer  function 
(i)V/3t  )  =  (  >V/P<i>  )  (  VI  /Tf  ) ,  where  (iV/3,  )  is  evaluated  for  a  dc  SQUID  with 

a  0  6a  6 

loop  inductance  biased  at  applied  flux  .  Similarly,  the  voltaqe  noise 
6  6 

spectral  density  is  S  ,  and  the  current  noise  spectral  density  through  the 

v  9  2 

inductance  Lp  is  S^  =  M  SJ/(Li  +  Lp),  where  Sy  and  Sj  are  evaluated  for  loop 
inductance  Le  and  applied  flux  4  .  The  correlation  factor  for  the  equivalent 
SQUID,  =  (S  S,  -  sf,)/S^,  is  identical  to  the  correlation  factor  for  the 

p  V  J  V  J  V  J 

isolated  SQUID  with  inductance  Lp  and  applied  flux 
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B.  DC  SQUID  with  Lossy  Input _ 

A  magnetometer  system  consisting  of  a  dc  SQUID  coupled  inductively  to 
an  input  circuit  containing  a  resistive  element  is  shown  in  Fig. 20(b).  A 
normal  resistance  is  in  series  with  the  pick-up  coil  inductance  Lp  and 
coupling  coil  inductance  Li-  The  effects  of  stray  capacitance  within  the  input  circ 
and  between  the  input  circuit  and  the  SQUID  are  not  included  in  the  analysis. 

The  introduction  of  loss  into  the  input  citcuit  destroys  the  quantization  of 
flux  within  that  loop.  As  a  result,  the  behavior  of  the  lossy  magnetometer 
differs  significantly  from  the  behvaior  of  magnetometers  with  superconducting 
input  circuits. 

The  equations  of  motion  for  the  system  are  as  follows.  The  equations 
describing  the  development  of  the  junction  phase  drops  (Eqn.  26-27)  and  the 
expression  for  the  voltage  across  the  SQUID  (Eqn.  25)  are  identical  to  the 
isolated  SQUID  equations.  The  equation  describing  fluxoid  quantization  in 
the  SQUID  loop  is  similar  to  the  superconducting  magnetometer  relationship. 


D  =  LJ  +  MH  +  f'g  (  39) 

Note  that  we  have  introduced  a  bias  flux  into  the  SQUID  loop.  In  the 
isolated  SQUID  operated  in  a  flux  locked  loop,  the  bias  flux  usually  consists 
of  both  dc  and  ac  components.  The  analog  of  Eqn.  33  is  qiven  by  the  voltage 
drop  around  the  input  circuit, 


(Lj+Lp)  dH/dt  +  MdJ/dt  +  R^  H+  d^/dt  +  VNJ  =  0  (40  ) 

where  V^j  is  a  random  voltage  generated  by  the  Johnson  noise  in  the  resistor 
Rj-  Eliminating  the  input  circuit  current  H  between  Eqn.  39  and  Eqn.  40 
yields 


d(D_LeJ)-  =  .  (D-LJ)  +  !B  +  d_h 

(Tt  i  i  dt 


(Li+L  )(dt  wl  ’ 


(41  ) 
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where  a  =  (L.  +  L  )/R^  and  the  screened  inductance  is  Lg  (Eqn.  36).  Thus 
the  equations  determininq  the  output  voltage  V  as  a  function  of  the  applied 
flux  >,•'  are  Eqn.  26,  Eqn.  27  and  Eqn.  41. 

d 

The  basic  equations  for  the  lossy  magnetometer  are  independent  of  the 
dc  value  of  the  externally  applied  flux  1  .  As  a  result,  the  system  cannot 
be  described  by  an  equivalent  device  with  forward  transfer  function  . 

d 

Instead,  we  shall  compute  the  response  of  the  system  to  a  time  dependent 
applied  flux  :  (t). 

d 

As  an  example  of  the  technique  to  be  used,  consider  the  following  case. 

Let  the  input  circuit  parameters  and  coupling  efficiency  be  such  that  the 

2 

screened  loop  inductance  Lg  L  (a  <<  1).  Consider  the  noise  free  case, 

=  0,  with  bias  flux  d :  g/d t  =  0.  Eqn.  41  reduces  to 

~  (D-LJ)  =  -  +  F (t )  +  :B/.  (42  ) 

where  F(t)  -  -  M (d :  /d t ) / ( l  ■  *•  L  ).  This  expression  can  be  integrated  for  a 

dip 

single  Fourier  component  of  F(t)  of  the  form  F^  sin  t 

D  -  JL  =  C  exp  (-t/j)  +  An  sin  (>nt  -  0n)  +  ^  (43) 


where 


An  =  Td+.aV)"4  Fn  (44) 

and 

0p  =  tan~^  (hiT ) .  (45  ) 

The  constant  C  is  determined  by  the  initial  conditions.  Note  that  a  sudden 
chanae  in  the  dc  bias  flux  l>g  introduces  a  transient  with  C  =  .V‘.g.  A  sudden 
change  in  the  amplitude  of  the  ac  applied  flux  at  time  t  introduces  a 
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transient  with  C  =  AAp  sin  (..•  tQ-  ")• 

Eqn.  43  is  the  analog  of  the  fluxoid  quantization  expression  for  the 
isolated  dc  SQUID  given  by  Eqn.  28.  The  form  of  the  two  equations  are  simi¬ 
lar.  The  lossy  input  circuit  has  the  effect  of  introducing  a  transient  re¬ 
sponse  time  into  the  quantization  condition.  A  sudden  change,  in  the 

dc  bias  flux  through  the  SQUID  loop  introduces  a  transient  M  exp  (-t/). 

In  addition,  an  ac  flux,  ,  applied  to  the  pick-up  coil  is  equivalent  to 
an  ac  flux  through  the  isolated  SQUID  loop  with  an  amplitude  and  phase  shift 
determined  by  the  input  circuit  parameters,  signal  frequency  .,  and  response 
time  i. 

The  well  known  solutions  for  the  isolated  dc  SQUID  may  now  be  exploited 

to  determine  the  effective  ac  forward  transfer  function  and  the  flux  resolution 

of  the  lossy  magnetometer  system.  Suppose  that  the  applied  flux  is  oscillating 

sinusoidally  at  frequency  .  •-<  «■,.  We  use  a  quasistatic  approximation 

to  determine  the  amplitude  of  the  output  voltage,  ,  at  the  signal  frequency  . 

in  terms  of  the  amplitude  of  the  applied  flux  I  and  the  isolated  SQUID 

dO 

forward  transfer  function,  :>V/  > ;■  , 

a 


v 


o 


a M(3V/3*  ) 

a 


R.  (l+i.2!2)* 


$ 

ao 


(46) 


At  frequencies  w  <<  1/x,  the  output  signal  is  reduced.  The  resolution  is 
determined  by  comparing  the  flux  noise  to  the  effective  signal  flux  in  the 
isolated  SQUID.  The  mean  squared  flux  noise  in  bandwidth  B  is  Af^  =  S^B, 
where  S,  is  evaluated  at  the  SQUID  bias  flux  (small  signal  limit).  The 
mean  squared  signal  flux  / 5  through  the  isolated  SQUID  is  related  to  the 

mean  squared  applied  flux  A!,  via 

a 


nl.l 
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(47) 


The  minimum  observable  applied  flux  with  1:1  signal-to-noise  ratio  (A>fs  =  At^) 
is  optimized  for  input  circuit  parameters  =  Lp. 


In  that  case 


f 


1 


Note  that  the  flux  resolution  deteriorates  at  frequencies  ,  ■  1/  like 

l/.";2. 

For  comparison,  the  minimum  mean  square  signal  flux  derived  from  the  CTG 

2 

model  for  an  input  coil  in  the  limit  R.  =  0,  a  «  1  i s 


;i 


(49) 


The  optimal  input  inductance  in  this  limit  is  L°pt  =  L  .  Thus  the  CTG  model 

^  P 

adequately  describes  the  behavior  of  the  lossy  magnetometer  in  the  weak 
coupling  limit  at  signal  freguencies  ..  -  l/i. 

Equation  48  was  derived  in  the  weak  coupling  limit  for  the  case  V^j  =  0. 
This  approximation  is  valid  provided  that  the  noise  rounding  of  the  magneto¬ 
meter  system  characteristics  are  dominated  by  the  Johnson  noise  sources 
associated  with  the  shunt  resistances  in  the  SQUID  rather  than  by  the  Johnson 
noise  in  the  input  circuit  resistance.  We  estimate  the  relative  importance 
of  the  two  noise  sources  as  follows.  The  input  noise  appears  in  Eqn.  41  as 
a  random  forcing  function  F  =  -  MVnj/((_.  +  Lp).  Since  Eqn.  41  is  linear,  F 
may  be  decomposed  into  a  Fourier  series  and  the  components  integrated 
following  the  method  used  for  the  applied  flux.  The  spectral  density  of  the 
resultant  flux  noise  in  the  integrated  equation  is 
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(50) 


where  Syj  is  the  voltage  power  spectral  density  of  the  Johnson  noise  generated 
by  the  input  resistance.  For  comparison,  the  independent  Johnson  noise 
sources  associated  with  the  junction  shunts  may  be  transformed  into  an 
equivalent  bias  current  noise  source  and  an  independent  flux  noise  sourceT' 

The  spectral  density  of  the  intrinsic  flux  noise  source  is 
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where  S^j  is  the  voltage  power  spectral  density  of  the  Johnson  noise  in  a 
shunt  of  resistance  R.  We  limit  the  discussion  to  system  parameters  such 


that  iij  1  / 1 .  In  this  case,  the  shunt  flux  noise  dominates  the  input 
equivalent  flux  noise  at  frequencies  m  in  the  optimized  system,  if 

SVI  «  (2/a2)(Lp/L)(<.l/R)2  SNI  (52) 

In  the  thermal  limit,  this  constraint  is  equivalent  to  the  condition 
RiVLp  <<  RT/L,  where  the  input  and  shunt  temperatures  are  T.  and  T,  and 
Wjl/R  s  i.  in  this  limit,  the  SQUID  characteristics  determined  by  the  total 
noise  spectral  density  near  the  Josephson  frequency  are  independent  of  V^j , 
and  the  input  Johnson  noise  source  may  be  approximated  as  a  low  frequency 
(quasi-static)  noise  source. 


V.  ANALYSIS  OF  DOUBLE  LOOP  MAGNETOMETER  RESOLUTION 

In  this  section  we  compute  the  magnetic  field  resolution  of  double 
loop  magnetometer  systems.  The  double  loop  SQUID  may  be  used  directly  as 
a  magnetometer,  or  coupled  to  an  input  circuit.  In  this  first  case,  the 
flux  resolution  of  the  device  in  the  limit  ;j  <<  a-j  is  approximately  equal 
to  that  of  a  conventional  SQUID  with  the  same  total  loop  inductance.  The 
double  loop  configuration  is  useful  in  those  cases  where  a  shortinq  capacitance 
across  the  SQUID  loop  is  unavoidable.  An  example  is  a  thin  film  gradiometer 
in  the  form  of  a  figure  eight  with  both  junctions  located  in  one  of  the 
loops.  In  the  second  case,  the  magnetometer  resolution  depends  on  both 
the  SQUID  and  input  circuit  parameters.  As  an  example,  we  discuss  the  reso¬ 
lution  of  a  double  loop  SQUID  coupled  inductively  to  a  superconducting  and 
normal  input  circuit  in  detail. 

A.  Double  Loop  SQUID  with  Superconducting  Input  Circuit 

The  magnetic  field  resolution  of  the  double  loop  SQUID  can  be  enhanced 
by  coupling  the  device  to  a  superconducting  input  circuit  (Fig.  21  (a)). 

The  circulating  currents  in  the  SQUID  oscillating  at  the  Josephson  frequency 
and  its  harmonics  generate  flux  in  the  input  circuit.  Screening  currents 
in  the  input  circuit  act  back  on  the  SQUID,  altering  the  behavior  of  the 
device.  As  a  result,  the  noise  spectral  densities  and  forward  transfer 
function  of  the  device  are  a  function  of  theinput  circuit  parameters. 

It  is  useful  to  model  the  double  loop  SQUID  plus  input  circuit  as  a 
single  superconducting  magnetometer  with  forward  transfer  function  and  noise 
spectral  densities  referred  to  the  pick-up  coil  inductance,  Lp  (See  Fig. 

21  (a)).  The  equivalent  magnetometer  can  then  be  used  to  replace  a  conventional 
dc  SQUID  in  magnetometer  and  linear  amplifier  applications  (see  Ref.  25). 

The  pick-up  coil  inductance  Lp  corresponds  to  the  conventional  SQUID  loop 
inductance  L.  The  forward  transfer  function  is  3V/T;  ,  where  t  is  the 

3  3  n 

external  flux  applied  to  the  pick-up  coil.  The  flux  noise  S,  =  S  /(">7/?<:  )  , 

'♦V  3 

and  the  noise  current  S„  through  L  correspond  to  the  conventional  flux 

*  *  r 

noise  and  circulating  current  noise  through  the  SQUID  loop  inductance  L. 


m 


The  equations  describing  the  double  loop  SQUID  plus  input  circuit  for 

the  superconducting  input  circuit  (Fig.  21(a))  are  as  follows.  The  total  flux 

applied  to  the  input  circuit  is  the  external  quasistatic  flux  t  plus  the 

2  °  2 

induced  flux,  MK,  generated  by  the  double  loop  SQUID,  where  M  =  a  L^L.. 

For  convenience,  we  shall  assume  that  the  nutual  inductance  between  L.  and 
L  is  negligible.  Since  flux  is  quantized  in  the  input  circuit,  the  screening 
current  H  satisfies 


H(Li  +  L  )  +  KM  +  *a  =  n:0  (53) 

In  addition,  the  fluxoid  quantization  condition  (Eqn.  23)  in  the  SQUID  becomes 
D  =  LJ  +  LSK  +  MH  (54) 

Thus,  Eqn.  23  is  replaced  by 


D  =  LJ  +  LeK  -  <e  (55) 

where 

r  =  L./(L.  +  Lp),  (56) 

't’e  =  ar  (Ls/Li  <fa  (57) 

and 

Le  =  Ls  (1  -  o 2r)  (58) 


The  equations  describing  the  development  of  the  phase  drops  ^  and  ^  (Eqn.  20- 
22)  are  unchanged.  Eqn.  24  becomes 


L  4  ♦  .  If 

e  dt7  s  dt 


(59) 


B 


Thus,  the  basic  equations  for  the  double  loop  SQUID  plus  input  circuit,  Eqns. 

20-22,  are  identical  to  the  original  equations  with  effective  coupling  loop 

2 

inductance  L  =  L  (1  -  «  r)  and  applied  flux  >:■  . 

e  s  e 

The  forward  transfer  function  which  relates  the  change  in  output  voltage 
AV  to  the  change  in  externally  applied  flux  AT  is  'W/"> ;■  =  (aV/at  )(jt  /%;  ) 

a  _  d  GG3 

and  the  voltage  noise  at  the  output  is  S° ,  where  V  and  Sy  are  computed  for  a 

double  loop  coupling  inductance  L  and  applied  flux  ■:  .  The  effective  flux 

e  o  2  e 

noise  through  the  pick-up  coil  is  S  =  S,c  L./-<  r  L  .  For  an  effective  total 

t't  1  s 

inductance  a  1,  ;>e  o,  <t>Q/4, 

Stp  ~  4kBT  0-  +  Ls(l  -  <t2r )12  Li./a2r2LsR. 


The  double  loop  magnetometer  plus  input  coil  can  now  be  represented  by  an 
idealized  magnetometer  with  input  coil  inductance  L  and  effective  flux  noise 
S.p. 

^  2 

The  flux  noise,  S  ,  is  minimized  in  the  perfect  coupling  limit,  .  =1. 

’  r 

In  that  case. 


S<*ip  '  (4kBT/R)  tLLi  +  LLp  +  LsLp^/LiLs-  (60) 

The  flux  noise  is  minimized  for  fixed  L  as  a  function  of  L.  as  follows: 

P  i 

LiPt  =  Lp  d  +  Ls/D  (61) 

S°f  3  (16  kBT/R)(l  +  L/Ls)LLp.  (62) 

We  are  interested  in  the  limit  L  L.  In  that  case,  L°pt=  IL  /(Lc  +  2L)  L 
and  =  Lp(l_s/L).  The  optimal  flux  noise  in  this  case  is  limited  by  the 

SQUID  loop  inductance  L  and  input  pick-up  coil  inductance  L  .  The  inductance 
L  is  constrained  by  Sj  2p  _  1.  (Note  that  the  case  B  =  =  1  has  been  solved 

numerically  in  Section  III.D.) 

The  current  noise  power  spectral  density  through  the  pick-up  coil  is 
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2  2 

Sr  =  M  S^/(L-  +  L  )  ,  where  S^  is  the  current  noise  power  spectral  density 

for  a  double  loop  SQUID  with  screening  loop  inductance  L  and  applied  flux 

2  e 

lg-  In  the  limit  L$  '  ■  t.,  -«  =  1,  Sr  S ^[-/ ( LS ke/^- p )  •  The  correlation  factor 
for  the  conventional  SQUID  is  Cp  =  (SySj  -  Syj)/Syj-  Note  that  the  correla¬ 
tion  factor  for  the  current  H  and  the  output  voltage  V  in  the  double  loop 
input  circuit  combination  is  identical  to  the  correlation  factor  for  the  current 
K  and  V  of  the  equivalent  double  loop  SQUID.  Furthermore ,  the  correlation 
factor  for  the  double  loop  SQUID  with  t  =  1,  =  1  is  comparable  to  the  cor¬ 

relation  factor  of  a  conventional  dc  SQUID  with  |t_  =  2. 

It  is  instructive  to  compute  the  energy  factor  referred  to  the  pick-up 
coil,  Sep  =  Spp/2Lp.  In  the  limit  L$  •  L,  Spp  ~  4kBT(2L)/R.  The  energy 
factor  of  a  conventional  dc  SQUID  with  loop  inductance  L^.  is  Sp  4kpT/R. 

This  is  a  reassuring  result.  The  noise  temperature  T^  of  a  conventional  SQUID 
linear  amplifier  with  a  tuned  input  circuit  operated  at  frequency  is  pro¬ 
portional  to  the  energy  and  correlation  factors,  Tf,  ~  -CpSp/KB25  The  optimal 
source  resistance  scales  with  the  input  SQUID  loop  inductance  L.  If  we  replace 
the  conventional  SQUIO  in  the  linear  amplifier  circuit  with  the  double  loop  plus 
input  circuit  combination,  the  energy  factor  is  unchanged  for  |_T  =  2L.  In 
addition,  the  correlation  between  the  output  voltage  and  the  current  noise 
in  the  pick-up  coil  is  similar  tot  he  correlation  between  the  SQUID  output 
voltage  noise  and  SQUID  loop  circulating  current  noise.  As  a  result,  the  optimal 
noise  temperatures  of  the  two  devices  are  similar.  This  is  important  because 
the  noise  temperature  of  a  high  gain  linear  amplifier  is  limited  by  the  uncer- 
tainty  principle  to  T^  ^  Wkgln  .  This  limit  is  approached  for  values  of 
the  device  parameters  such  that  Kmj  '■->  kgT.  In  this  case,  the  thermal  noise 
sources  in  the  shunted  junctions  are  dominated  by  quantum  fluctuation  effects 
in  the  SQUID.  Minimum  noise  temperatures  computed  as  a  function  of  the  con¬ 
ventional  dc  SQUID  parameters  are  consistent  with  the  uncertainty  principle 
requirement^  Since  the  behavior  of  the  conventional  dc  SQUID  and  the  optimized 
double  loop  SQUID  plus  input  circuit  with  Lp  =  21  are  similar,  we  expect  a 
detailed  analysis  would  also  yield  Sp  fi ,  Cpp  Cp  in  anreement  with  the 
uncertainty  principle  limit.  However,  the  optimal  source  resistance  of  the 


I 


double  loop  amplifier  now  scales  with  the  pick-up  inductance  Lp.  This  value 
is  not  limited  by  the  uncertainty  principle. 

The  procedure  for  optimizing  the  magnetic  field  resolution  is  the  follow¬ 
ing.  The  flux  noise  for  the  double  loop  magnetometer  can  be  reduced  by 
lowering  the  junction  loop  inductance  L.  The  coupling  inductance  is  increased 
to  L$  s>  L  to  improve  the  coupling  properties  of  the  double  loop  SQUID.  The 
inductance  L  is  then  screened  by  the  superconducting  input  circuit.  The 
actual  input  inductance  for  the  combined  system  is  Lp  =  (L./Ls)L,  where  L.  is 

the  largest  coil  that  can  be  efficiently  coupled  to  the  inductance  L  .  The 

2  s 

inductance  screening  technique  is  only  advantageous  in  the  limit  u  1.  For 
o'  ?  1,  L  '•'»  L,  the  effective  inductance  L  [L  +  2L  (1  -  a^)]/(2  -  ) . 

S  0  5 

If  the  effective  inductance  is  dominated  by  L  rather  than  by  L,  the  effective 

2  ^ 

flux  noise  deteriorates  like  (L$/a  L).  This  places  a  limit  on  the  values  of  L, 

L$  and  L. ,  and  determines  the  value  of  Lp.  The  magnetic  field  resolution  is 

then  a  function  of  the  geometry  of  L  ,  and  S^-p. 

One  of  the  properties  of  the  double  loop  SQUID  is  that  it  permits  the  use 

of  SOUID  oeometries  with  improved  coupling  properties.  For  example,  a  thin 
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film  planar  double  loop  SQUID  has  been  fabricated  by  Jaycox  and  Ketchen.  In 
addition,  a  toroidal  point  contact  SQUID  similar  to  a  double  loop  SQUID  has 
been  fabricated  by  Paik?^  In  both  cases,  the  coupling  properties  of  these 
devices  are  superior  to  the  conventional  planar  and  cylindrical  dc  SQUIDs.  As 
a  result,  the  double  loop  SQUID  is  a  more  promising  device  for  the  investigation 
of  inductance  screening  effects.  Note,  however,  that  the  analysis  has  ignored 
any  stray  capacitance  in  the  circuit.  The  effects  of  such  a  capacitance  must 
be  negligible  even  at  the  effective  Josephson  frequency. 


B.  Double  Loop  SQJID  with  Lossy  Input  Circuit 

The  analysis  of  the  lossy  double  loop  magnetometer  follows  the 

analysis  of  the  lossy  dc  SQUID  system  closely.  The  double  loop  SQUID 

is  coupled  inductively  to  the  input  circuit  as  shown  in  Fig.  21(b). 

The  mutual  inductance  between  the  input  and  SQUID  coupling  loop  is 
2  2 

M  =  a  L$.  The  mutual  inductance  between  the  input  circuit  and 
the  SQUID  loop  inductance,  L,  is  taken  to  be  negligible,  as  is  the 
effect  of  any  stray  capacitance  in  the  circuit.  The  pick-up  coil 
inductance  is  L„  and  the  loss  is  FL . 


The  voltage  drop  around  the  input  circuit  is 

<Li  tLp>$tM$+RiHtTr*vNi'°-  (63) 

where  the  flux  $  is  applied  to  the  pick-up  inductance  L  ,  H  is 
a  p 

the  current  flowing  in  the  inout  coil,  and  the  Johnson  noise  in  the 
loss  is  V^j.  Quantization  of  the  fluxoid  in  the  SQUID  yields 

D  =  LJ  +  LK  +  MH  +  4>0,  (64) 

S  D 

where  the  bias  flux  to  the  SQUID  loop  is  Eliminating  the  circulating 

current  H  in  the  input  circuit. 


iL(D-LeK-U)=- 


(D-LK-LJ)  *  dOR 
- 5 -  +  —  +-Hr 

T  j  T-  dt 


where  x.  =  (L.  +  L  )/R.  as  before,  and  the  screened  coupling  loop 
l  l  p  1  2  _  i 

inductance  Lg  =  t_s(l  -  a  r)  for  r  *  L.(L  +  L.)  .  The  voltage 

drop  around  the  coupling  loop  in  the  SQUID  yields  the  analog  of 

Equation  24, 


i: 


L 


The  junction  equations,  Eqn.  20-21,  and  the  output  voltage  relationship, 
Eqn.  22,  are  unchanged. 

Following  the  analysis  for  the  conventional  lossy  magnetometer, 

2 

we  consider  the  weak  coupling  limit  (ct  «  1,  Lg  ~  L  ).  In  this 
case  Eqn.  65  has  the  form 


3F  0  '  L,  K  -  LJ>  ■  - 


(0  -  L  K  -  LJ) 


T  . 
1 


+  F(t), 


(67) 


where  the  forcing  function  F(t)  is  a  function  of  the  input  and  applied 
flux,  and  the  input  Johnson  noise  voltage.  Integrating  Eqn.  67  for  a 


single  Fourier  component  of  F(t)  of  the  form  F^  sin  t, 

D  -  L  K  -  LJ  =  C  exp(-t/i.)  +  A  sin(w  t  -  6  ), 
s  r  l  n  n  n 


where 


and 


A  =  x.(l  +  w2  F  , 

n  v  n  'i'  n’ 


'n  '  tan''  Ti>' 


(68) 

(69) 

(70) 


Combining  Eqn.  68  and  Eqn.  66  for  constant  bias  flux  4^, 

2 

M  ^4  =  — -  exp( - 1/ 1  . )  -  u)2  A  sin(w  t  -  0  ). 

dt  4 

Eqn.  71  can  be  used  to  eliminate  the  current  in  the  input  coil,  H, 
from  Eqn.  66  .  For  times  t  t.,  the  transient  response  can  be 
neglected.  We  express  Eqn.  66  in  terms  of  an  effective  circulating 
current  Ke> 

d2K  dK  K  ,  c!V 

L  _ e  +  d  _e.  +  _e  =  R  _d J.  +  jL  +  _J. 

Ls  Ks  dt  C  s  dt  C  dt  ’ 


(71) 


(72) 
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JMiu, 


'  ***  ** 


where 


(73) 


where 

K 

=  K  - 

■  K  sin(w 

- 

O') 

» 

(73) 
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o  n 

n 
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ro 

C  A  [(1  - 

2 

w 

L 

C  )2  +  (u>  R  C  )2) 

(74  ) 

0 

n 

s  n 

n 

s 

s  K  n  s  s 

and 

l 

0 

n 

=  °n 

+  tan  ^  [m 

n 

R 

s 

Cs 

<>-“n  Ls 

(76  ) 

2  i 

In  the  low  frequency  limit,  oin  «  1/L  Cs>  u>n  <<  1/R$  C^,  9n  -  0n 
and  K0  =  K  -  Cs  An  sin(w|(t  -  0p). 

The  fluxoid  quantization  condition,  Eqn.  68  ,  can  also  be 
expressed  as  a  function  of  the  effective  circulating  current  Kg  for 
times  t  »  i.  , 


D  =  L  K  +  LJ  +  4>  ,  (76  ) 

s  e  e 

where  the  effective  flux 

°e  =  An  ~  °n}  +  Ls  Ko  sin((V  ’  0n>‘  (77  1 

Eqn.  76  and  Eqn.  77  now  have  the  same  form  as  the  corresponding 
equations  for  the  isolated  double  loop  SQUID.  As  a  result,  the  complete 
set  of  equations  for  the  lossy  double  loop  magnetometer ,  Eqns.  20,  21, 

22,  76  ,  and  77  ,  can  be  solved  using  the  same  techniques  for  the 
isolated  driven  double  loop  SQUID.  For  example,  consider  a  lossy 
double  loop  magnetometer  system  with  parameters  satisfying 
ft .  t .  »  ft  f  t .  >> 1  and  R  C  ft  «  1  (see  Sec.  Ill  B) .  The  effective 
flux  noise  components  of  4>e  generated  by  the  input  resistance  at 
frequencies  >  ft  will  produce  noise  rounding  of  the  double  loop 
SQUID  characteristics.  In  this  limit,  e'  =  0  .  The  effective  flux 
noise  spectral  density  near  resonance  is  S^j  ~  Syj/R^  >  where  Syj 
is  the  voltage  power  spectral  density  of  the  Johnson  noise  generated 
by  the  loss  R. .  The  input  and  SQUID  parameters  must  be  chosen  such 
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that  the  smearing  of  the  characteristics  produced  by  S  ,j  is  negligible 
(see  Sec.  III.D).  This  places  an  upper  limit  on  the  input  resistance 
and  temperature 

The  effective  flux  noise  at  frequencies  ta  ^  SI,  is 
p  n  J 

S  I  ■-  4S WT/w  •  The  equivalent  flux  noise  aenerated  within  the  SQUID  is 
t  1  2  *  *  n  o 

S,s  =  L  S^j/2R  ,  where  S^j  is  the  voltage  spectral  density  of  the  Johnson 
noise  qenerated  by  the  junction  shunt  (see  Sec.  IV. B).  Thus,  the  noise 
rounding  of  the  device  characteristics  can  be  neglected  provided 


SVI  8  (ftj  L/R)2  SNI  ,  (78) 

where  the  factor  (q,  L/R)  1  for  SQUID  loop  g  =  1.  Note  that  if 
the  input  noise  cannot  be  neglected,  the  resultant  equations  cannot  be 
well  approximated  by  those  of  an  isolated  symmetric  double  loop  SQUID. 

We  assume  that  the  effects  of  the  input  Johnson  noise  can  be 
neglected.  In  that  case,  for  signal  frequency  w  <  Qt,  the  applied 
flux  $>  through  the  magnetometer  system  may  be  approximated  by  a 

a 

quasistatic  effective  flux  <?  through  an  isolated  double  loop  SQUID 
(see  Sec.  IV. B).  In  this  limit,  the  contribution  of  the  screening 
current,  K  -  «e>  to  the  effective  applied  flux  can  be  neglected.  As 
a  result,  the  analysis  of  the  resolution  follows  the  pattern  of  the 
conventional  lossy  dc  SQUID  magnetometer  (Sect.  IV. B).  The  amplitude 
of  the  output  voltage  generated  by  an  applied  flux  t  =  t>ao  sin  „st 
is 


V 


Ws  M  ( 9  V  /  3  4>a ) 


ao 


Ri(l 


2x1/2 
f  i ) 


sin(u  t 


■6)  > 


(79) 


where  u  =  tan’*(u  v-),  and  (jV/at  )  is  the  quasistatic  foreward 
transfer  function  for  the  isolated  double  loop  SQUID.  The  minimum 
observable  mean-squared  flux  At,  in  bandwidth  B  for  unity  signal- 

Cl 

to-noise  ratio  is 


At 

a 


+-s 

I,- 


B 
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for  optimal  input  inductance  L^°*3+  =  L  •  *n  exp>ession,  S 
flux  noise  spectral  density  of  the  isolated  double  loop  SQUID. 


VI.  SUMMARY 

An  analysis  of  double  loop  SQUID  magnetometer  systems  has  been  performed. 

The  device  consists  of  a  conventional  dc  SQUID  with  loop  inductance  in  sene  , 
with  a  resonant  coupling  circuit.  A  lumped  circuit  element  model  for  the 
device  has  been  developed.  The  equations  of  motion  for  the  device  in  the  noise 
free  case  are  seen  to  be  closely  related  to  the  conventional  dc  SQUID  equations 
for  devices  with  coupling  loop  resonant  frequency  well  below  the  intrinsic 
Josephson  frequency.  The  I-V  characteristi cs  are  found  to  be  hysteretic,  with 
critical  current  determined  by  the  total  loop  inductance.  In  addition,  possibly 
chaotic  behavior  is  observed  for  devices  with  resonant  frequency  on  the  order 
of  the  Jospehson  frequency. 

Noise  rounded  characteristics  were  generated  by  direct  numerical  simula¬ 
tion.  Optimal  flux  resolution  is  obtained  at  bias  currents  on  the  order  of 
the  total  loop  critical  current.  The  device  characteristics  at  that  point 
are  well  approximated  by  conventional  dc  SQUID  characteristics  for  a  device 
with  loop  inductance  equal  to  the  double  loop  SQUID  junction  loop  inductance. 
However,  the  hysteresis  in  the  double  loop  characteristics  limits  the  minimum 
flux  resolution  to  a  value  which  is  well  approximated  by  the  resolution  of  a 
conventional  device  with  the  same  total  loop  inductance. 

The  main  advantage  of  the  double  loop  SQUID  appears  to  lie  in  the  enhanced 
coupling  properties  which  the  double  loop  geometry  permits.  With  that  in  mind, 
a  model  for  SQUID  magnetometer  systems  which  explicitly  includes  the  effects 
of  the  input  circuit  on  the  SQUID  characteristics  was  developed.  The  model 
was  used  to  calculate  the  resolution  of  superconducting  and  lossy  magnetometer 
system.  Inductance  screening  effects  obtainable  for  input  circuits  coupled 
inductively  to  the  SQUID  were  found  to  enhance  the  magnetic  field  resolution 
in  the  good  coupling  limit.  Deterioration  of  the  flux  resolution  of  lossy 
systems  was  observed  at  signal  frequencies  below  the  inverse  of  the  input 
circuit  decay  time.  The  effects  of  stray  capacitance  within  the  input  circuit 
and  between  the  input  circuit  and  the  SQUID  were  not  considered  in  this  effort. 
These  effects  could  be  substantial  if  resonant  behavior  in  the  input  circuit 
is  observed  at  the  couplina  loop  resonant  frequency  or  at  the  Josephson  frequency. 
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Figure  5.  Noise  Rounded  1-V  Characteristics  for  a  DC  SQUID 
as  a  Function  of  Applied  Flux 


Figure  6. 


Low  Frequency  Circulating  Noise  Spectral  Density 
S°/2r  Versus  Bias  Current  I/IQ  as  a  function 
of  ft  for  a  DC  SQUID. 
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Low  F requency  Correlation  Spectral  Density  S  ,/2r 
Versus  Bias  Current  I / 1  as  a  function  of 
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Figure  12.  Average  Voltage,  •:  v>,  Versus  Applied  Flux 
(>j..  for  a  DC  SQUID  with  5=1. 
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Noise  Free  I-V  Characteristics  for  a  Double  Loop 
SQUID  with  8*1,  8^  3  50  and  $>a  =  0.25. 
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Figure  15.  Noise  Free  I-V  Characteristics  for  a  Double  Loop 
SQUID  with  p,  =  1,  flL  *  10,  and  <j,a  =  0.25  as 
a  function  of  the  tank  circuit  frequency  Wf. 
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Figure  16.  Plot  of  Voltage  v,  and  Screening  Currents  k 
and  j  as  a  Function  of  Time,  o  for  a  Double 
Loop  SQUID  with  3  =  1,  p.  -  10.  Pq  0.3537, 


Figure  18.  Noise  Rounded  I-V  Characteristics,  (b)  Forward 
Transfer  Function,  (c)-(g)  Low  Frequency  Noise 
Power  Spectral  Densities,  and  (h)  Dimensionless 
Flux  Noise  as  a  Function  of  Bias  Current  for  a 
j  Double  Loop  SQUID  with  $  =  1,  eL  =  1,  3^  =  100, 

p  =  0.01  and  r  =  0.01 
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Figure  19. 


Noise  Rounded  I-V  Characteristics,  (b)  Forward 
Transfer  Function,  (c)-(g)  Low  Frequency  Noise 
Power  Spectral  Densities,  and  (h)  Dimensionless 
Flux  Noise  as  a  Function  of  Bias  Current  for  a 
Double  Loop  SQUID  with  s  =  1,  0^  =  10,  0^  =  10, 
p  =  0.1  and  r  =  0.001. 
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